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ABSTRACT 


Research  under  this  Grant  has  developed  new  theoretical  and  computa¬ 
tional  methods  for  studying  the  properties  of  materials.  In  particular,  it  has 
addressed  the  problems  of  using  ab  initio  computations  to  predict  (a)  chemi¬ 
cal  bonding  in  systems  composed  of  heavy  elements,  (b)  relativistic  effects  in 
molecules,  and  (c)  core-electron  spectra,  such  as  obtained  by  XPS  experiments, 
in  heavy  atoms  and  molecules.  Practical  computational  techniques  have  been 
formulated  and  tested  which  allow  not  only  improved  understanding  of  the 
theory  of  materials,  but  also  routine  calculations  of  energies  and  structures  of 
a  variety  of  atomic  and  molecular  systems. 
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AFOSft  77-9148 

DEVELOPMENT  &  PRACTICAL  APPLICATION 
OF  THE  THEORY  OF  MATERIALS 


I.  Introduction 

Any  Air  Force  which  does  not  keep  its  doctrines  ahead  of  its  equipment  and 
if*  visions  far  into  the  future  can  only  delude  the  nation  into  a  false  sense  of 
security. 

General  H.  H.  Arnold 


The  research  carried  out  under  Grant  AFOSR-77-3145  has  examined 
some  of  the  most  basic  principles  governing  the  ways  in  which  materials  are 
composed  in  teems  of  atoms  and  molecules.  Our  first  objective  was  to  pursue 
the  concept  of  core  electrons,  which  are  of  comparatively  high  energy  and  do 
not  participate  directly  in  bonding  between  atoms,  as  opposed  to  the  lower- 
energy  valence  elctrons  that  do  tend  to  form  bonds.  This  distinction  allows 
extremely  efficient  ab  initio  computational  methods  for  the  calculation  of 
molecular  structure,  such  as  the  NOCOR  procedure  we  developed  during  the 
previous  Grant  period.  For  example,  we  have  shown  that  an  increase  in  com¬ 
putational  efficiency  of  roughly  an  order  of  magnitude  for  molecules  composed 
of  elements  from  successive  rows  of  the  Periodic  Table  may  result. 

Our  original  method  was  largely  limited  to  nonrelativistic  single¬ 
configuration  calculations  using  simple  basis  sets  to  describe  the  electron  dis¬ 
tributions.  This  was  found  to  be  a  surprisingly  accurate  way  of  describing  rela¬ 
tively  large  molecular  systems  composed  of  moderately  heavy  elements,  such  as 
phosphorus,  silicon  and  germanium.  However  during  the  present  Grant  period 
we  have  extended  both  the  theory  and  computer  programs  to  make  the  proce¬ 
dure  much  more  general.  Thus  it  is  now  applicable  to  any  set  of  basis  func¬ 
tions,  either  single-  or  multiconfiguration  calculations,  species  with  unpaired 
electrons,  and  has  the  facility  to  treat  relativistic  effects,  which  predominate  in 
many  properties  of  heavy  elements  such  as  lead  and  uranium.  Thus  we  have 
now  developed  a  computational  tool  for  predicting  a  variety  of  spectroscopic 
and  chemical  properties. 
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A  number  of  discoveries  in  related  fields  were  an  outgrowth  of  this  effort. 
From  a  theoretical  viewpoint,  by  far  the  most  significant  was  the  develop¬ 
ment  of  the  first  completely  relativistic  theory  of  molecular  structure,  along 
with  the  associated  computer  programs.  For  this  we  have  employed  the  self- 
consistent- field  approach  using  the  Dirac  (relativistic)  equation  of  electron  mo¬ 
tion.  Implementation  of  the  Dirac  equation  involved  considerable  computa¬ 
tional  complexity,  since  all  orbitals  are  four-component  complex  functions  in¬ 
stead  of  one-component  real  functions  as  in  the  nonrelativistic  approximation. 
Although  our  program  is  presently  at  a  level  only  analogous  to  the  restricted 
Hartree-Fock  method,  we  believe  it  may  well  provide  a  cornerstone  of  sub¬ 
sequent  work  in  other  laboratories  on  the  theoretical  chemistry  of  heavy  ele¬ 
ments.  Relativistic  effects  are  now  known  to  be  quite  significant  in  determining 
the  chemical  properties  of  elements  from  the  fifth  and  lower  rows  of  the  Periodic 
Table.  For  example,  the  (relativistic)  spin-orbit  energies  of  lead  compounds  are 
often  greater  than  their  bond  energies.  Therefore  an  a b  initio  technique  for 
predicting  such  effects  is  of  central  importance  to  work  in  this  field. 

Another  outgrowth  of  our  study  of  core  electron  properties  and  relativis¬ 
tic  corrections  was  a  procedure  for  calculating  high-energy  ionization  spectra, 
such  as  X-ray  photoelectron  spectra  (XPS),  of  atoms.  This  type  of  spectroscopy 
measures  the  minimal  energy  required  to  remove  an  electron  from  a  core  orbi¬ 
tal.  We  found  that  employing  previously-developed  Dirac-Fock  techniques  for 
atoms,  a  striking  agreement  with  the  available  experimental  results  could  be 
attained.  For  example,  differences  between  theory  and  experiment  are  often 
only  3%  of  the  lowest  errors  previously  reported  -  an  error  which  approaches 
the  experimental  linewidth.  We  have  now  found  this  agreement  to  be  consis¬ 
tent  for  core  electrons  with  experimentally  accessible  energies  throughout  the 
Periodic  Table. 

These  are  the  three  principal  areas  to  which  we  have  made  contributions 
during  the  tenure  of  this  Grant.  Others,  which  have  evolved  from  our  im¬ 
plementation  of  general-purpose  a b  initio  programs,  will  be  discussed  below. 
In  the  following  sections  we  will  very  briefly  review  the  work  published  under 
this  Grant.  We  then  summarize  work  that  has  been  completed  but  not  yet 
published,  and  our  conclusions.  In  the  Appendix  are  some  selected  papers  that 
represent  the  central  theoretical  and  conceptual  achievements  under  this  Grant. 
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II.  Summary  of  Work  Published  Under  This  Grant 


A.  Theoretical  Advances 

As  we  applied  our  original  nocor  approach  to  systems  composed  of  in¬ 
creasingly  unusual  and  heavy  elements,  it  soon  became  apparent  that  we  would 
have  to  confront  the  problem  of  relativistic  effects  -  i.e.  those  embodied  in 
the  special  theory  of  relativity  and  manifested  due  to  the  high  velocity  of  core 
electrons  in  the  region  of  space  near  a  nucleus  of  high  nuclear  charge.  These 
types  of  terms  affect  essentially  every  type  of  molecular  property,  including 
bond  lengths  and  heats  of  formation.  At  the  time  we  began  this  work,  it  was 
not  known  how  large  these  effects  would  be.  However  we  and  others  have  sub¬ 
sequently  shown  that  they  can  be  significant  in  systems  of  practical  interest. 
Our  first  paper  on  this  subject1  showed  how  our  nocor  procedure  has  a  precise 
and  theoretically  rigorous  analog  in  the  Dirac-Fock  formalism.  We  proved  the 
interesting  result  that  all  electrostatic  terms  coupling  core  and  valence  orbitals 
in  the  effective  core  potential  vanish  except  for  one  (small)  term.  We  developed 
a  Dirac-Fock  procedure  for  calculating  valence-only  molecular  properties  under 
the  assumption  that  only  core  electrons  move  at  relativistic  velocities.  Finally, 
we  showed  how  a  practical  computational  procedure  using  this  approach  could 
be  developed,  and  reported  sample  calculations  on  the  spectroscopic  properties 
of  the  PbO  molecule.  In  this  compound,  since  there  is  only  one  heavy  element, 
the  relativistic  shift  is  fairly  modest:  about  0.01  A  in  the  bond  length  and  10 
cm”1  in  the  vibrational  frequency.  A  copy  of  this  paper  is  included  below  as 
Appendix  A. 

Our  next  paper2  in  this  series  dealt  with  the  formation  of  effective  poten¬ 
tials  at  a  basic  theoretical  level.  We  demonstrated  how  our  effective  potential 
procedure  could  be  extended  to  a  somewhat  higher  level  of  accuracy  by  employ¬ 
ing  distinct  core  potentials  for  each  valence  orbital.  This  in  turn  allows  the 
construction  of  rigorous  Phillips-Kleinman  pseudopotentials  for  systems  with 
unpaired  electrons  and  ,  more  importantly,  for  the  case  of  configuration  mixing. 
For  the  latter  we  have  developed  an  effective- potential  multiconfiguration  SCF 
program.  This  procedure  allows  the  calculation  of  all  molecular  properties  that 
depend  strongly  on  the  electron  correlation  or  reflect  a  breakdown  in  the  single¬ 
configuration  model,  such  as  bond  energies  and  the  detailed  shapes  of  potential 
energy  surfaces.  With  some  modifications  and  restrictions  it  also  allows  the 
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calculation  of  the  energy  separation  of  electronic  states.  As  an  example  we 
showed  in  some  detail  how  the  method  works  for  the  Li2  molecule,  which  is 
small  enough  that  a  parallel  all-electron  calculation  is  feasible.  The  results  of 
the  two  calculations  in  terms  of  valence  orbitals  and  dissociation  energies  were 
found  to  be  in  excellent  agreement.  A  copy  of  this  paper  is  included  below  in 
Appendix  B. 

B.  nocor  Studies  of  Representative  Systems 

At  the  same  time  we  were  extending  the  theoretical  basis  of  our  computa¬ 
tional  procedure,  we  continued  to  apply  our  existing  programs  to  molecular 
systems  that  exhibit  unusual  and/or  theoretically  important  types  of  chemical 
bonding.  For  example  we  have  studied  the  structure  of  the  polyiodide  ions3, 
Ig“  and  Ig~.  These  have  a  polymeric  structure  that  is  not  manifested  by  similar 
chemical  elements.  In  a  crystal  lattice,  they  are  known  to  have  an  unusual  zig¬ 
zag  shape  with  sets  of  three  1  atoms  being  colinear.  We  calculated  the  detailed 
electronic  and  geometrical  structure  of  these  ions.  In  addition  we  performed  an 
extensive  study  of  the  way  in  which  the  presence  and  symmetry  of  the  crystal 
lattice  can  influence  the  ionic  geometry.  We  found  the  remarkable  result  that 
in  the  I~  compounds  the  charge  alternates  between  even-  and  odd-numbered 
atoms  in  a  regular  fashion.  In  addition  the  bending  force  constants  also  exhibit 
a  very  marked  alternation,  with  every  second  bond  angle  being  stiff,  and  alter¬ 
nate  angles  having  a  force  constant  close  to  zero  allowing  the  angle  to  vary 
freely.  Thus  in  vacuo  or  in  a  symmetric  crystal  lattice,  these  chains  would  be 
linear.  However  in  a  lattice  that  is  even  slightly  asymmetric,  the  chains  are  bent 
in  the  experimentally  observed  pattern.  This  is  apparently  a  type  a  chemical 
behavior  that  had  not  been  observed  before. 

The  next  molecule  of  this  series  that  we  studied  was  the  organometallic 
compound  cyclopentadienylthallium(I)4.  The  structure  of  this  molecule  is  a 
“half-sandwich”  with  the  bond  from  the  metal  to  the  cyclopentadiene  ring  being 
symmetric  about  the  center  of  the  ring.  This  type  of  bonding  in  a  heavy  ele¬ 
ment  such  as  thallium  had  not  been  theoretically  studied  before.  We  performed 
a  detailed  analysis  of  its  electronic  structure  to  determine  the  exact  nature  of 
the  bonding.  We  found  that  about  one  half  of  the  bonding  is  electrostatic  in¬ 
teraction  between  a  T1+  and  C5H5"  ion.  The  covalent  bonding  was  found  to  be 
primarily  of  ir  symmetry,  as  might  be  expected.  Further,  although  the  neutral 
T1  atom  has  an  electronic  configuration  s2p,  which  is  not  favorable  for  such 


Final  Report  for  AFOSR-77-314S 


7 


bonding,  in  the  molecule  almost  a  full  electron  charge  (0.69  e)  is  promoted  to 
contracted  outer  d  orbitals.  This  strongly  favors  x  bonding.  Specifically  most 
of  the  metal-ring  bonding  results  from  a  valence  orbital  of  ei  symmetry  in  which 
the  thallium  atom  has  pd2  hybridization.  This  hybridization  is  significant  in 
that  it  is  a  “directed  x  bond”,  that  is,  of  an  e  symmetry  but  strongly  polarized 
in  one  direction.  From  this  observation  it  was  possible  to  explain  the  structure 
of  not  only  this  compound,  but  many  chemically  related  ones.  Thus  the  metal¬ 
ring  bond  is  weaker  in  CsHsIn,  which  lacks  available  valence  d  orbitals.  The 
elements  Au  and  Hg  lack  the  valence  p  electron  and  hence  cannot  form  the  type 
of  bond  we  observed  in  the  T1  compound.  This  explains  why  these  elements 
form  a  rather  than  x  complexes.  Thus  our  calculation  was  able  to  correlate 
much  of  the  experimental  data  on  these  systems. 

The  next  three  papers  of  this  series  were  done  in  collaboration  with  Mr. 
Michael  Rothman,  a  graduate  student  from  the  University  of  Michigan,  Ann 
Arbor,  who  came  to  this  site  to  learn  the  NOCOR  method.  After  returning  to 
that  campus  he,  along  with  Professors  Lawrence  Bartell  and  Lawrence  Lohr, 
began  a  productive  research  program  independent  from  ours,  using  our  com¬ 
puter  programs  to  study  unusual  inorganic  systems.  We  describe  here  only  the 
work  in  which  we  were  directly  involved. 

The  first  study  we  proposed  in  this  series  dealt  with  the  Te2-*-  ion5.  This 
is  an  unusual  system  in  several  respects:  it  is  an  inorganic  aromatic  ring  struc¬ 
ture  that  satisfies  the  4n+2  rule,  and  has  a  Hiickel  bond  order  of  1.5.  It  was  one 
of  the  first  Te  compounds  to  be  discovered,  is  extremely  stable  in  solution,  and 
displays  a  remarkable  bright  red  color.  We  undertook  a  detailed  study  of  its 
geometric  and  electronic  structure,  with  and  without  the  presence  of  oppositely- 
charged  counterions.  We  obtained  a  bond  length  in  vacuo  of  2.6108A,  and 
2.6739A  with  counterions.  The  latter  compares  well  with  the  experimental  bond 
length  in  crystals  of  2.6633  -  2.6738A.  The  molecular  structure  was  found  to  be 
square  and  planar,  as  expected.  We  were  also  able  to  calculate  the  electronic 
transition  that  accounts  for  the  color.  This  was  the  first  practical  application 
of  our  new  effective- potential  mcscf  program.  Although  it  had  been  the  object 
of  much  debate  previously,  we  were  able  to  show  that  this  transition  is  x  -♦  x* 
and  we  were  also  able  to  assign  another  experimentally  observed  transition  as 
n  -*  x*.  Finally  we  calculated  all  of  the  (five)  vibration  frequencies.  Four  of 
these  had  been  assigned  experimentally.  We  were  able  to  prove  that  two  of  the 
four  vibrational  bands  had  previously  been  misassigned  due  to  being  masked 
by  stronger  lines  in  the  experimental  spectrum. 
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The  final  two  papers  of  this  series  concerned  the  properties  of  the  hyper- 
valent  xenon  fluorides.  The  first  paper6  studied  the  structures  and  detailed 
intramolecular  potential  energy  surface  of  XeF2  and  XeF4.  Although  these 
compounds  are  of  some  interest  in  their  own  right,  we  also  felt  it  to  be  desirable 
to  perfect  the  basis  set  and  other  aspects  of  the  computation  before  attacking 
the  more  highly  fluorinated  compounds.  A  minimal  basis  set  was  developed 
that  is  considerably  more  accurate  for  this  type  of  calculation  than  the  usual 
atom-optimized  set.  The  calculation  reproduced  the  correct  symmetries  of  both 
compounds.  The  bond  lengths  that  resulted  were  somewhat  too  long,  but  gave 
the  correct  difference  between  the  two  molecules.  Also,  although  the  stretch¬ 
ing  force  constants  were  slightly  too  large,  the  interaction  force  constants  were 
found  to  be  in  good  agreement  with  experiment,  as  was  the  anharmonicity. 
This  was  the  first  time  that  such  a  detailed  analysis  of  the  vibrational  potential 
of  such  a  relatively  large  molecule  had  been  carried  out.  We  felt  it  reflected 
sufficient  accuracy  to  warrant  the  study  of  the  two  additional  xenon  fluorides, 
whose  structures  are  of  considerably  more  theoretical  and  experimental  inter¬ 
est. 

These  two  systems,  which  were  the  subject  of  the  subsequent  paper7,  are 
XeFj^  and  XeF6-  Both  of  these  are  of  considerable  theoretical  and  experimen¬ 
tal  interest,  since  they  have  long  been  known  to  exhibit  a  number  of  puzzeling 
thermodynamic  and  spectroscopic  properties.  The  structures  of  both  species 
are  governed  by  a  stereochemically-active  lone  pair  of  electrons.  We  undertook 
a  very  detailed  study  of  the  structures  of  both  species  and  also  their  force  fields, 
including  quadratic,  cubic  and  quartic  terms.  The  structure  we  obtained  for 
XcF+  was  found  to  be  in  excellent  agreement  with  experiment.  For  example 
the  angle  between  bonds  to  the  axial  and  equitorial  F  atoms  was  calculated 
to  be  80.8°,  in  comparison  to  the  experimental  value  of  79.0°.  However  the 
major  question  has  long  been  the  structure  of  XeF6.  It  is  known,  from  electron- 
diffraction  data,  to  be  nonrigid  and  possess  a  Csv  symmetry  that  is  slightly 
deformed  from  an  octahedron.  Our  calculations  bore  this  out,  showing  the 
structure  of  lowest  energy  to  be  of  this  symmetry.  Plotting  the  electron  density 
in  the  molecular  orbitals,  we  were  able  to  show  graphically  the  presence  and 
position  of  the  lone  electron  pair.  This  protrudes  through  one  of  the  triangular 
faces  if  the  structure  is  constrained  to  be  octahedral,  and  hence  causes  its  dis¬ 
tortion.  The  final  structure  we  obtained  is  in  good  agreement  with  electron 
diffraction  results:  the  angle  between  the  Csv  axis  and  an  F  atom  in  the  opened 
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face  was  79.4°,  which  may  be  compared  to  the  experimental  value  of  67.5°. 
The  angle  between  the  C3W  axis  and  an  F  in  the  opposite  face  was  47.4°  in 
comparison  the  the  experimental  value  52.3°.  However  the  barrier  to  hindered 
pseudorotation,  that  is  the  energy  to  move  the  electron  pair  between  different 
faces  of  the  molecule,  was  somewhat  too  high  to  explain  the  rapid  interconver¬ 
sion  of  conformations  that  is  observed  experimentally  at  room  temperature. 
This  reflects  the  fact,  shown  by  the  above  distortion  angles,  that  our  calculated 
structure  was  slightly  too  distorted.  Nonetheless,  our  calculations  ended  a  long 
controversy  by  confirming  the  structure  indicated  by  spectroscopic  data  and 
demonstrating  its  electronic  origins.  Finally,  in  comparing  XeF6  and  XeF^“ 
with  the  other  two  xenon  compounds,  it  was  found  that  the  trends  in  their 
respective  orbital  energies  and  equilibrium  bond  lengths  were  well  described  by 
our  calculations. 

C.  Ab  Initio  Vibrational  Spectra 

As  an  outgrowth  of  our  work  on  unusual  types  of  bonding  and  their 
confirmation  by  computed  spectroscopic  properties,  we  calculated  the  infrared 
and  Raman  spectra  of  the  thiirene  molecule8,  C2H2S.  This  project  was  carried 
out  in  collaboration  with  Professors  Schaad  and  Hess  of  this  department. 
Thiirene  is  a  ring  compound,  which  due  to  its  strongly  antiaromatic  electronic 
structure  is  unstable  and  may  be  studied  only  by  its  spectra  in  matrices  at  low 
temperature.  These  types  of  spectra  invariably  contain  lines  from  many  other 
species.  Thus  it  is  extremely  useful  to  calculate  those  parts  of  the  spectrum 
which  correspond  to  a  particular  specie*  in  order  to  confirm  its  structure.  For 
this  work  we  employed  an  all-electron  SCF  calculation  using  a  moderately  large 
(4-31 G)  basis.  We  computed  all  the  vibrational  bands  of  both  C2H2S  and 
C2D2S,  including  both  their  frequency  and,  in  the  case  of  infrared  bands,  their 
relative  intensities.  Our  computed  structure  exhibited  an  unusually  long  C- 
S  bond  and  a  fairly  short  (1.2509A)  C=C  bond,  as  might  be  expected  from 
the  system’s  antiaromatic  electronic  structure.  We  were  able  to  reproduce  all 
the  lines  of  the  experimental  spectra  in  each  species  except  for  a  pair  of  lines 
in  each,  which  we  showed  had  been  experimentally  misassigned.  The  spectros- 
copists  subsequently  reevaluated  the  spectra  and  brought  the  assignments  into 
accord  with  our  calculations.  The  final  result  of  this  work  was  to  confirm  the 
postulated  structure  of  this  species. 
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D.  Mechanism  of  the  H2  Exchange  Reaction 

During  our  wor*  with  MCSCF  calculations  on  small  molecules,  we  became 
interested  in  the  H2  +  D2  exchange  reaction®.  This  system  is  important  from 
a  theoretical  viewpoint  as  the  prototype  for  all  four-center  reactions  and  also 
of  the  energy  surfaces  of  reacting  diatomic  molecules  in  general.  However  for 
a  number  of  years  this  reaction  has  completely  defied  all  theoretical  explana¬ 
tion,  since  the  experimental  evidence  indicates  a  four-center  transition  state 
with  an  activiation  energy  much  less  than  any  that  has  been  calculated  from 
first  principles.  In  collaboration  with  Michael  Rothman  and  Lawrence  Lohr 
of  the  University  of  Michigan,  Ann  Arbor,  we  began  a  series  of  calculations 
of  interacting  H2  molecules  using  a  large  basis  set  and  a  set  of  20  optimized 
configurations  in  an  MCSCF  calculation.  We  developed  a  technique  for  studying 
this  type  of  system  by  which  the  centers  of  mass  of  the  interacting  molecules  are 
gradually  brought  together,  all  other  degrees  of  freedom  being  simultaneously 
optimized  to  provide  the  minimum  energy.  Along  this  path  we  then  searched 
for  the  transition  state,  defined  as  the  point  at  which  the  energy  could  decrease 
along  one  of  these  other  degrees  of  freedom,  thus  dissociating  into  products.  In 
this  way  we  located  an  energy  which  is  a  lower  bound  to  the  transition  state 
energy.  Our  value,  193  kJ  mol-1  is  34  kJ  mol-1  above  the  experimental  ac¬ 
tivation  energy.  Hence  we  were  able  to  demonstrate  that,  at  least  to  this  level  of 
accuracy,  the  experimental  measurements  are  inconsistent  with  a  bimolecular 
mechanism. 

E.  Computation  of  Core  Energy  Levels  in  Atoms  and  Small  Molecules 

As  an  outgrowth  of  our  work  on  core-valence  partitioning  and  relativis¬ 
tic  effects,  we  began  a  very  productive  collaboration  with  M.  S.  Banna 
of  this  department,  who  works  in  the  area  of  X-ray  photoelectron  spectra 
(XPS)  of  atoms  and  small  molecules.  This  type  of  spectroscopy  employs  a 
monochromatic  beam  of  X-rays  to  ionize  core  electrons,  measuring  the  mini¬ 
mum  energy  required  for  the  ionization  to  take  place.  The  XPS  spectra  of  free 
atoms  pose,  except  for  the  inert  gases,  a  difficult  experimental  problem,  but  one 
that  is  central  to  understanding  the  interaction  of  high-energy  radiation  with 
matter.  It  is  also  useful  in  studying  the  nature  of  bonding  in  solid  materials. 
We  realized  that  our  computer  programs  might  be  very  useful  in  predicting 
and  interpreting  the  XPS  spectra  in  two  cases:  small  molecules  composed  of 
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light  atoms,  for  which  relativistic  effects  are  small;  and  heavy  atoms,  for  which 
the  relativistic  effects  are  comparatively  well- understood.  Surprisingly  good 
accuracy  was  found  for  the  core  ionization  energies  we  calculated  in  both  types 
of  systems.  This  has  been  the  basis  of  three  papers  so  far. 

In  the  first  paper10,  Dr.  Banna  reported  for  the  first  time  the  X-ray 
spectrum  of  an  alkali  halide  in  the  gas  phase,  Csl.  To  check  and  interpret  his 
results,  and  also  to  test  our  computational  procedure,  we  carried  out  relativistic 
(Dirac-Fock)  calculations  of  the  observed  3d5/2  ionization  in  Cs,  I,  and  their 
singly-charged  ions.  Using  these  results  in  conjunction  with  the  experimental 
energies,  it  was  possible  to  analyze  the  nature  of  the  bonding  of  the  Csl  molecule 
in  terms  of  covalent  and  ionic  contributions.  We  also  studied  the  ionization 
process  in  terms  of  atomic  and  extraatomic  relaxation  energies,  which  also 
reflect  the  charge  distribution  of  the  neutral  molecule.  One  of  the  most  useful 
results  of  this  work  was  the  discovery  that  the  relativistic  atomic  structure 
program  gave  very  precise  estimates  of  core  binding  energies  for  experimentally 
known  transitions.  They  were  also  appreciably  more  accurate  than  the  best 
theoretical  values  previously  available.  For  example,  the  best  previous  value  for 
the  ionization  of  the  Cs  atom  was  749.7eV.  We  calculated  732.1,  while  the  ex¬ 
perimental  value  is  731.8.  Likewise  for  the  I  atom,  the  best  previous  value  was 
643. OeV.  We  calculated  627.9,  while  the  experiment  gives  628.1.  The  linewidth 
of  the  X-ray  line  is  typically  O.leV.  Therefore  we  concluded  that  our  method 
is  sufficiently  accurate  that  it  may  be  used  to  calculate  ionization  in  that  great 
majority  of  atomic  systems  that  can  not  be  studied  experimentally,  and  these 
results  may  be  reliably  used  to  interpret  experimental  results  for  molecules. 

Our  second  paper11  dealt  with  the  use  of  nonrelativistic  molecular  SCF 
programs  to  study  the  ionization  of  small  molecules  composed  of  light  elements. 
We  compared  the  compounds  ozone,  O3,  and  oxygen  difluoride,  OF2.  These 
have  long  been  an  enigma  in  the  understanding  of  the  electronic  structures 
of  these  types  of  molecules,  because  the  ionization  of  the  central  O  atom  of 
ozone  is  of  very  high  energy  compared  to  other  oxygen  compounds,  indicating 
a  large  positive  charge  on  that  atom.  It  is  higher  than  the  oxygen  of  OF2,  even 
though  F  is  a  much  more  electronegative  element.  We  calculated  all  four  of 
the  core  ionization  energies,  using  a  large  (4s3p)  gaussian  basis  set  and  finding 
the  energies  of  the  open-shell  hole  states  using  the  BISONMC  MCSCF  program. 
First  we  showed  that  the  calculation  gives  good  agreement  with  experiment, 
except  for  the  central  O  of  ozone.  However  all  the  experimentally  observed 


1 2  Final  Report  for  AFOSR-77-31 45 

trends  were  reproduced.  The  error  for  the  central  O  of  ozone  was  1.68eV.  For 
comparison,  the  error  for  the  terminal  O  was  0.04eV.  In  OF2  the  error  for  the 
O  atom  was  0.23  and  for  the  F  atoms,  0.44.  In  terms  of  relative  accuracy, 
these  final  three  values  are  in  error  by  approximately  0.05%  .  We  calculated 
the  central  oxygen  of  ozone  to  possess  an  ionization  2.45eV  higher  than  in  OF2; 
the  experimental  value  is  l.OeV.  This  difference  is  probably  due  to  our  neglect 
of  d  orbitals  in  the  basis  set  for  ozone,  since  these  are  known  to  to  affect  the 
central  atom  most  strongly.  We  then  carried  out  a  population  analysis  of  the 
electronic  structures  of  both  compounds.  This  showed  that  although  the  F 
atom  withdraws  considerable  charge  from  the  O  atoms  in  OF2,  it  does  so  only 
through  the  o  orbitals.  The  transfer  is  small  through  the  ir  orbitals  even  though 
the  7r  overlap  is  large.  The  transfer  of  charge  is  small  in  the  o  orbitals  of  ozone, 
but  very  large  and  in  the  opposite  direction  in  the  ir  orbitals.  Approximately 
0.2  electron  is  shifted  from  the  central  oxygen  in  this  way.  Thus  the  XPS  spectra 
essentially  reflect  the  f  bonding  present  in  ozone  but  not  in  OF2. 

In  our  third  paper12,  we  carried  out  a  series  of  comparisons  of  the  atomic 
core  electron  ionization  energies  of  several  atoms  and  atomic  ions.  The  choice 
of  the  particular  elements  and  orbitals  chosen  was  dictated  by  the  availability 
of  experimental  data.  We  compared  Hartree-Fock,  Dirac-Slater,  and  our  Dirac- 
Fock  approaches.  In  this  case  we  also  included  J-dependent  terms,  thus  ob¬ 
taining  all  four  possible  transitions  that  would  correspond  to  a  spin-orbit  pair 
or  a  single  Hartree-Fock  orbital.  Our  results  were  found  to  be  more  accurate 
than  any  that  have  been  previously  reported,  and  this  agreement  is  seen  to  be 
largely  independent  of  atomic  number  or  valence  electron  configuration.  This 
paper  is  included  below  as  Appendix  C. 
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III.  Summary  of  New  Theories  and  Computer  Programs 

A.  Comparison  of  Effective  Potential  Methods  of  Varying  Levels  of  Accuracy 

As  noted  at  the  beginning  of  this  report,  our  original  NOCOR  method  for 
carrying  out  valence-only  molecular  calculations  was  straightforward  to  apply 
to  any  system  and  gave  quite  accurate  results,  but  was  limited  to  nonrelativis- 
tic  single-configuration  calculations.  Now  we  have  removed  these  restrictions. 
So  it  is  of  interest  to  examine  how  their  presence  or  absence  affects  molecular 
properties.  We  therefore  undertook  a  study  of  the  spectroscopic  properties  of  a 
quite  heavy  diatomic  compound,  Pb2.  To  obtain  reasonable  accuracy,  we  first 
needed  to  generate  a  basis  set,  since  the  standard  tables  do  not  extend  this  far 
down  the  Periodic  Table.  We  therefore  carried  out  a  number  of  atomic  calcula¬ 
tions  to  optimize  a  complete  double-zeta  basis  for  the  Pb  atom.  The  valence 
part  (6s,6p)  is  then  suitable  for  forming  molecular  orbitals.  We  then  performed 
nonrelativistic  scf  and  six-configuration  MCSCF  calculations  to  determine  the 
total  valence  energy  Ew,  equilibrium  bond  length  Re,  first  vibrational  frequency 
we  and  anharmonicity  t uexe.  We  then  calculated  these  same  properties  by  both 
the  SCF  and  mcscf  methods  but  employing  relativistic  (Dirac-Fock)  core  orbi¬ 
tals.  The  results  are  shown  and  compared  with  the  available  experimental  data 
in  Table  I.  The  experimental  Re  is  not  known,  but  is  thought  to  be  about  2.9  A. 
From  this  table  it  may  be  seen  that  the  relativistic  valence  energies  are  much 
lower  than  the  nonrelativistic  ones,  reflecting  the  strong  contraction  of  the  s 
and  p  core  orbitals  in  the  former  case.  SCF  bond  lengths  are  always  shorter 
and  vibrational  frequencies  higher  than  in  the  more  accurate  MCSCF  results, 
since  SCF  wavefunctions  must  dissociate  to  spurious  excited  electronic  states 
at  infinite  bond  length.  Also  relativistic  calculations  are  always  more  tightly 
bound,  with  shorter  bond  lengths  and  higher  vibrational  frequencies  due  to  the 
contraction  of  the  core  orbitals.  The  effect  of  correlation  energy,  present  in 
the  MCSCF  results,  is  about  as  large  as  those  due  to  relativity.  But  the  final 
results,  embodying  both  effects,  are  the  most  accurate.  Any  calculation  that 
employed  one  but  not  the  other  of  these  effects  would  have  implied  an  incorrect 
conclusion  as  to  which  terms  are  important  in  determining  the  nature  of  the 
bonding.  This  work  is  currently  being  prepared  for  publication. 

B.  Full  Dirac-Fock  Theory  of  Molecules 
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In  order  to  incorporate  relativistic  effects  rigorously  into  effective  poten¬ 
tial  calculations,  it  is  first  necessary  to  first  develop  a  method  for  determining 
the  exect  four-component  (Dirac-Fock)  orbitals  for  molecules.  This  had  not 
heretofore  been  achieved,  due  to  computational  complexity.  In  our  previous 
work  we  had  employed  only  relativistic  core  orbitals,  thus  avoiding  the  calcula¬ 
tion  of  multicenter  terms  with  these  types  of  orbitals.  However  it  is  clearly 
desirable  to  treat  the  valence  on  this  level  also,  since  it  would  allow  the  use 
of  j-j  coupling  and  hence  such  important  properties  as  spin-orbit  splitting  in 
molecules.  We  have  now,  in  collaboration  with  Prof.  S.  N.  Datta  who  was 
a  postdoctoral  associate  with  our  group,  developed  a  complete  new  computer 
program  which  achieves  this  result.  This  program  solves  the  Dirac-Fock  equa¬ 
tions  for  both  atoms  and  molecules  using  a  basis  set  of  gaussian  functions, 
which  are  appropriate  for  computing  multi-center  integrals. 

It  is  interesting  that  in  the  Dirac  theory,  the  total  energy  of  an  atom  or 
molecule  is  unstable  in  a  variational  calculation  with  respect  to  the  formation 
of  anti  particles,  in  this  case  a  positrons,  with  the  decrease  in  energy  of  2mc*. 
This  is  very  similar  to  the  variational  collapse  that  occurs  in  valence-only  cal¬ 
culations  when  the  pseudopotential  is  not  present.  Datta  was  able  to  prove 
that  our  effective  potential  (NOCOR)  formalism  could  be  applied  directly  to  the 
case  of  variational  collapse  in  the  Dirac  equation13.  This  is  the  only  method 
derived  thus  far  which  permits  Dirac-Fock  calculations  on  molecules  without 
the  variational  collapse  to  -2rac2,  and  has  been  included  in  our  program. 

In  determining  appropriate  basis  sets,  it  was  soon  discovered  that  entirely 
different  procedures  were  required  from  the  nonrelativistic  case.  For  example 
the  imaginary  terms  proved  much  more  difficult  to  expand  than  the  real  ones. 
Also,  due  to  the  presence  of  two  coupled  radial  components,  increasing  the 
number  of  basis  functions  did  not  always  give  a  more  precise  result.  This 
problem  was  solved  by  constraining  the  basis  always  exactly  to  reproduce  the 
correct  one-electron  terms  in  the  relativistic  energy.  It  was  found  that  ten-term 
gaussian  basis  sets  gave  a  good  estimate  of  the  relativistic  correction  to  the 
energy  and  were  sufficiently  compact  for  molecular  calculations.  For  example, 
using  this  procedure  to  generate  a  basis  for  the  Be  atom,  we  found  a  total  energy 
of -14.574078  hartrees.  The  analogous  nonrelativistic  value  is  -14.572579.  The 
exact  relativistic  value  found  using  numerical  methods  is  -14.575189.  Since  the 
known  relativistic  energy  correction  is  .00217,  it  is  seen  that  our  simple  ten- 
term  basis  has  given  69%  of  this  effect.  We  have  calculated  the  Be2  molecule  in 
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the  tame  way  and  found  essentially  the  same  qualitative  result.  Of  course  Be2 
is  too  light  to  have  a  measurable  shift  in  its  binding  energy,  i.e.  the  difference 
between  the  molecule  and  its  constituate  atoms. 

Our  first  paper  describing  this  work  is  nearly  finished.  We  hope  in  the 
near  future  to  extend  our  calculations  to  heavier  systems.  The  program  in 
essentially  its  present  form  may  be  used  in  effective- potential  calculations  as 
well,  where  the  effects  we  have  studied  will  be  much  more  pronounced. 

C.  Symmetries  of  Transition  States  in  Isoenergetic  Reactions 

Finally,  one  additional  finding  developed  from  our  work  on  the  diatomic 
exchange  reaction  described  earlier.  This  was  a  study  of  the  symmetries  of 
these  reactions  in  terms  of  their  nuclear  geometries.  In  the  absence  of  a  detailed 
knowledge  of  the  precise  geometry  of  the  transition  state  of  a  reaction,  one  may 
study  the  activation  energy  of  a  reaction  in  two  ways:  either  by  bringing  the 
centers  of  mass  of  two  interacting  molcules  together  along  the  minimal-energy 
pathway,  or  by  chooeing  the  minimal  symmetry  the  transition  state  must  have, 
and  then  minimising  the  energy  of  this  structure.  In  the  paper  described  above, 
we  took  the  first  approach.  But  in  the  course  of  this  work,  we  also  developed 
the  second,  which  may  be  at  least  as  useful. 

This  approach  is  restricted  to  the  special  case  of  isoenergetic  reactions,  i.e. 
those  in  which  products  and  reactants  have  the  same  total  energy.  In  this  case 
we  were  able  to  prove  that  the  transition  state  must  have  at  least  one  element 
of  symmetry.  Further  we  proposed  the  concept  of  'symmetrically  equivalent 
groups",  and  using  these  worked  out  the  symmetries  of  transition  states  for 
any  three-  and  four-center  reaction.  In  particular  we  proved  that  for  any  four- 
center  mechanism  involving  diatomic  molecules,  the  transition  state  must  be 
of  C2«  symmetry  with  a  perpendicular  symmetry  plane  passing  through  two, 
and  only  two,  of  the  atoms.  Thus  in  determining  the  transition  state  energy 
of  reacting  diatomics  by  ab  initio  computations,  one  need  only  optimise  the 
energy  of  a  structure  of  this  symmetry. 

Since  this  work  is  not  closely  related  to  any  of  the  others,  it  is  included 
below  as  Appendix  0. 
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The  Phillips  Klemman  pseudopolcntial  and  the  opeuloi  equivalence  techniques  for  calculating  effective  potentials  ere 
applied  directly  to  the  Dirac  -Hartrce-Foclc  calculation  of  molecular  properties.  The  results  based  on  both  relativistic  and 
nonrelativistic  core  functions  are  compared  for  the  PbO  molecule. 


I.  Introduction 

The  recent  development  of  effective-potential  theory  has  allowed  accurate  calculations  to  be  performed  on  the 
stiucture  and  properties  of  molecules  composed  of  elements  from  all  regions  of  the  periodic  table.  It  has  been  sug¬ 
gested  (I)  that  relativistjc  effecis  might  be  significant  in  compounds  of  fifth-row  and  heavier  elements,  such  as  Pb. 
Since  the  advantages  of  effective-potential  theory  in  quantum-chemical  calculations  should  be  most  pronounced 
for  compounds  of  such  heavy  elements,  it  would  be  especially  valuable  to  be  able  to  include  these  relativistic  con¬ 
tributions  in  the  core  potential. 

As  part  of  our  ongoing  investigations  of  molecular  quantum  calculations  involving  heavy  atoms,  we  conceived 
about  two  years  ago  the  approach  described  in  this  paper.  This  treatment  is  based  on  the  assumption  that  only  the 
core  orbitals  3rc  inherently  relativistic,  and  hence  relativistic  effects  on  the  valence  orbitals  ate  indirect,  through 
core-valence  electron-electron  interactions.  In  this  paper  we  develop  a  valence-electron  theory  based  on  the 
Dirac -Fock  hamiltonian.  specifically  defining  the  relativistic  Phillips-Kleinman  operator,  as  well  as  calculating 
the  corresponding  local  Coulomb  and  exchange  potential  functions.  Finally  we  apply  our  theory  io  some  proper¬ 
ties  of  the  PbO  molecule  and  compare  our  results  with  the  analogous  nonrelativistic  results,  with  experiment  and 
with  a  nonrelativistic  all-electron  calculation  (2). 


2.  Valence-electron  treatment  of  the  Dirac- Hart ree-Fock  theory 

The  relativistic  hamiltonian  for  a  many-electron  system  may  be  written  as  (3) 

//°  =  L//D<ji)+  £  (‘> 

M 

where  is  the  nuclear -nuclear  repulsion  term,  if  any,  and  //p(jr)  is  the  standard  Dirac  one-election  hamiltonian 
u  the  field  of  nuclear  potential  V, 

nuclei 

*vtl.  =>  -  £  Zp/r„ 

The  terms  in  eq.  (2)  retain  their  usual  meaning  (3|.  In  this  paper  we  will  not  consider  Breit  interactions  or  any 


(2) 
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other  higher-order  relativistic  corrections,  so  that  the  total  energy  may  be  written  as  the  expectation  value  of//0 


«<//«> £(./,* -K^+V". 


/,  -  W.ltfiW*  -<^0)^(2)lrfJ,|^(l)^(2)>.  -  <*#(l)^(2)|rfJ,|*#(2)^(l)>. 

If  we  assume  that  the  orbitals  of  the  system  may  be  divided  into  strongly-orlhogona!  core  and  valence  sets  (4), 
minimization  of  <//°>  is  equivalent  to  the  minimization  of  £°  ♦  where  Km,  is  the  screened  nuclear-nuclear 
repulaion  and  £°  it  given  by 

(«> 

where 

N  v  Ni,Ni  aue’ei 

//?«£  H'd(h)  ♦  J  £  £  r-  >  //J>(p)  -  //D( H)  ♦  £  5  (■'cp  -  Kcp) 

m  0+»  '  tr* 

and  Af£  are  the  conventional  direct  and  exchange  operators  involving  the/th  core  orbital  of  thepth  atom,  all 

core  orbitals  being  hela  constant. 

The  atomic  orbitals  that  minimize  the  value  of  an  atomic  </#°>  are  characterized  in  genera)  by  the  quantum 
numbers 

«  *  n,  ♦  |*|,  /•-(«*♦$).  «■-/.-/♦  I.-./  *nd  *■*! 
as  discussed  elsewhere  |3)  and  can  be  written  in  terms  of  the  two-component  angular  functions  ♦)  « 

_ _ i  \  -« 


e  n  _1  \ 

r.*,4)m'  w 


where  the  radial  functions  Pmt0)  and  QHK(')  are  conventionally  called  the  “large"  and  “small"  components.  The 
orthonormality  conditions  of  the  atomic  orbitals  are  given  by  (after  using  the  symmetry  properties  of  xNM) 

an 

/  *(£./*  ♦  Q.0*)  -  <are*«V4«..,s*  (6> 

0 

Next  we  coiMider  the  valence-electron  problem  for  an  atomic  case  In  terms  of  these  functions.  The  atomic  cote 
la  assumed  to  consist  of  a  set  of  relativtstlcslly  closed  groups  A,  each  with  the  electronic  population  q*  s  ty*  *  1. 
Following  the  treatment  of  previous  authors  (3]  one  can  show  that  for  atoms,  by  analogy  with  the  all-electron  case. 


*?-£/;♦!  ££  £  M*(«, »F*(«. fl  -  «*(o. »C*(o, «] , 

a* I  an*  k 

p[v«. 


/  dr(£#(c(«K?0/dr  -  «.0./r)  ♦  (c*  -  n#J  -Q.lc(dF./dr  ♦  %mPJr)  ♦  (rl  ♦  KX?J}. 

where  A*.  J*  and  are  constant  coefficients  a*  defined  by  previous  authors.  The  integrals  t\(a.  b)  and 

G*(a,  *)  are  gvan  by 
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Fk(a,  b)  =  J  l^(r)  +  Q* (rJlr-1  Yk(b,  b,r)dr , 

0 

m 

Gk(a,b)=  f  [Pa{r)Pb(r)  +  Qa(r)Qb(r)]r-'Yk(a,b-,r)<lr,  (8) 

0 

Yk{a,  b,r)  =  rf  +  Qa(s)Qb(s)\is  . 

0 

r<  and  r>  are  defined  as  the  lesser  and  greater  respectively  of  r  and  s.  The  variational  problem,  therefore,  corre¬ 
sponds  to  the  minimization  of  A®  as  given  in  eq.  (7)  subject  to  the  orthonormality  constraint  (6).  Thus  the  varia¬ 
tional  condition  is  of  the  general  form 

A£-v°  =  E  ea„  a[/  (/■„/•,  +  Qa  Qfi)ir  j .  (9) 

One  may  apply  a  unitary  transformation  to  the  radial  components  of  both  P„  and  Qa  such  that  the  lagrangjan  multi 
pliers  eag  are  brought  into  diagonal  form  ea(J  -*  €ft(J6nj4  The  re,ullin8  pseudoeigenvalue  equation  is 

written  in  the  matrix  form  as  a  0  a 


where  the  four  elements  of  the  F  matrix  are  given  by 

Fn  *(cJ-  f)+£[r-'y0(i<,j(;r)-|E  k,^(lk0,  (r)] 

t  EE  YkV,P.r)-B*{a,fi)r-*nk,(r)\  . 

0  k 

F{  2  ~  c(d/ir  -  K/r),  -  -c(d/dz  +  */f) , 

F 22  -  -(c*  +  n  +  E  [r->  E  O*  ,  (r)] 


+  E'E  lAH*,P)r-lYkQ},e,r)-B*(a,fi)r-'nk,(r) J . 

0  k 

Hie  primed  summations  are  over  all  functions  except  a.  The  operator  Slky  is  defined  by 

After  a  little  manipulation  one  obtains  the  familiar  form  of  the  valence-electronic  energy 
*v 

£?  -  E  (/;  ♦  £,)  , 

a»l 


(12) 
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where  Ca  =  eaa  can  be  considered  as  the  analogous  orbital  energy. 


3.  Relativistic  Phillips-Kleinman  pseudopotential  and  atomic  local  potential 


Since  eq.  (10)  is  of  t)ie  same  form  as  that  which  occurs  in  the  nonrelaiivistic  case,  wc  define  ihe  onc-electron 
form  of  the  atomic  Phillips-Kleinman  pscudopotential  [  5  J ,  after  integrating  over  the  angular  coordinates, 


=S  |(^)>(eo  -  e„  Q  ft)  |  X  angular  terms 

The  pseudoeigcnvalue  equation  (10)  now  reduces  to 


(13) 


(14) 


where  Pa  and  Qa  are  the  radial  components  of  the  ath  valence  pseudowavefunction. 

In  order  to  calculate  the  relativistic  analogue  of  the  atomic  local  potential  [6],  we  define  a  new  formal  Fock 
matrix  Fin  the  space  spanned  by  the  pseudowavcfunclions  as 


^U-EE  [Ak(a.(i)r-lYkU},t),r)  -  Bk(a,0)r-'Slkfi(r)]+(c2  -  F) , 

fi  K 

^IJ^II'  ^21  “  FH  • 

Fn  -  E’E  (J*(a,0)r-'yk(fi,  0,r)  -  Bk(a,0)r-'nkfi(r)}  -(c*  ♦  V) , 


(15) 


where 


rk(0.P;r) 


Qa(r) 


Pa(') 

Qa(r) 


m 

r  /  (r*/r*+,)[?fl2(r)t^(i))<lr 
0 


and 


>[?.<*)%(*>  ♦ 


Eq.  (14)  is  now  equivalent  to 


(F+Kn‘  +  R'«) 


(16) 


where  W®  is  an  atomic  local  potential  matrix,  and  +  IV®  =  Vt{t  is  the  effective  potential  matrix  fot  the  va¬ 
lence-electron  treatment.  Comparing  eqs.  (14)  and  (16)  one  can  write  from  the  operator  equivalence  condition  [6) 


Wa 


(17) 


whereby  IV®  can  be  explicitly  written  as 
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a  U/co  +  U/Va 

11.22  Wll . 22  +  *11,22  • 

h'ir.22-p«-  */ar-‘n^  (o]?0.ea/?..Ca  . 

0  x 

-Bk(°.P)'- 1 1 nkf(r)  -  nw(r))]?„  &/?„.&• 

Thus  after  determining  the  effective  potential  one  may  then  solve  the  pseudoeigenvalue  equation  (16)  for  an  atom 
and,  assuming  the  Wa  are  transferable,  for  a  molecule  [6]. 

4.  Nonrelativist ic  limit  for  the  valence  electrons 

Since  the  magnitude  of  the  relativistic  one-electron  effect  depends  on  the  shielded  nuclear  charge  for  any  orbital, 
it  is  reasonable  to  believe  that  the  valence-electronic  behaviour  is  approximately  nonrelativistic  except  through  the 
indirect  effect  of  interaction  with  the  core  electrons  for  whom  the  relativistic  effect  is  predominant.  Eq.  (10)  is 
equivalent  to  a  set  of  coupled  equations 

F\  lpa  *F\lQam  z*Pa>  F2lQ*  *  F22Pa  “  • 

Recognizing  ea  *  cJ  +  and  neglecting  the  terms  independent  of  c  in  the  second  coupled  equation  gives  the 
familiar  nonrelativistic  limit  as  v/c  -*  0,  where  u  is  the  electron  velocity.  A  contraction  of  the  modified  Fock  ma¬ 
trix  in  the  subspace  spanned  by  the  “large"  components  gives  the  pseudoeigenvalue  equation 

FPa  =  j  <dW)  +  Ka(K„  +  1)/2a*  ♦  £  ^  [/■-»  Y0(* ,  A ;  f)  -  i  £  ry-  kiar- 1  fi*  *  «] 


+  E'  ^  M*(a,/J)r-1  r)  -  5* (a,  ©r»  n^toljw  -  t'aPa(r)  . 


is  the  analogue  of  the  nonrelativistic  valence-orbital  energy,  where  the  core  has  been  described  relativistically. 
The  contracted  effective  potential  H'con  is  equal  to  IVj |  since  IVjj  has  been  neglected  in  comparison  to  c. 

A  further  simplification  arises  by  noting  that  at  c  -» •»,  Q  -*  0.  In  the  limit  C  =  0  we  will  denote  the  correspond¬ 
ing  quantities  with  a  prime.  In  that  case  the  orthogonality  condition  is  given  by 


With  this  choice  of  the  large  components  the  Phillips-  Kleinman  pseudopotential  is  strictly  deteimined  from  only 
the  large  components  of  the  core.  The  atomic  local  potential  is  also  considerably  simplified  and  is  defined  only  for 
the  large  component  of  a  pseudovalence-function.  This  is  then  identical  with  IV (  j  as  defined  by  eq  l  IS)  with  ihe 
necessary  modifications  of  the  radial  components,  Pa  -» Pa,  Qa  -*  0.  Thus  in  the  expansion  of  h'(  i  m  eq.  ( 1 8)  all 
the  Y  and  fl  integrals  will  be  independent  of  any  Q  components  in  this  case  except  the  integral  T0(y{ ,  sf ;  r),  in 
which  the  core  Q  components  still  contribute.  Therefore  one  may  generate  the  effective  potential  for  each  of  the 
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constituent  atoms  of  a  molecule  and  finally  solve  the  molecular  pseudoeigenvalue  problem  [6J.  The  molecular 
Pock  operator  so  obtained  is  defined  in  the  space  spanned  only  by  the  valence  pscudowavefunctions  corresponding 
to  the  radial  components?'  and  is  identical  in  form  with  its  nonrelativistic  analog 


S.  Sample  computations  and  discussion 

As  a  test  of  this  theory  we  have  carried  out  parallel  relativistic  and  nonrelativistic  calculations  for  some  proper¬ 
ties  of  PbO.  The  relativistic  atomic  functions  were  obtained  by  the  numerical  Dirac-Fock  procedure  [7 J .  The  va¬ 
lence  basis  functions  and  nonrelativistic  core  functions  were  taken  from  atomic  double-zeta  SCF  calculations 
The  valence  functions  were  defined  to  be  the  6s  and  6p  orbitals  of  Pb  and  the  2s  and  2p  orbitals  of  O.  The  core 
potentials  WCQ  and  the  pseudopolentiuls  were  found  by  the  method  described  above.  The  lv'v“  potential,  due  to 
valence -valence  interactions,  was  generated  from  the  nonrelativistic  atomic  orbitals  and  assumed  to  be  the  same 
in  both  cases.  The  Wa  functions  were  found  separately  for  each  member  of  the  basis  set.  For  the  O  atom,  only  non¬ 
relativistic  functions  were  employed. 

In  table  1  we  list  the  orbital  energies,  equilibrium  bond  lengths  and  vibrational  frequencies  obtained  from  these 
two  calculations.  We  also  compare  the  latter  two  properties,  for  which  the  SCF  approximation  usually  gives  accu¬ 
rate  results,  with  the  experimental  values  (9).  The  orbital  energies  are  seen  to  be  considerably  different  in  the  two 
calculations,  the  relativistic  values  being  lower,  although  the  bond  lengths  and  frequencies  are  nearly  the  same.  As 
might  be  expected  [  1  ] ,  the  bond  length  in  the  relativistic  case  is  slightly  shorter. 

In  interpreting  these  results  it  should  be  pointed  out  that  although  most  of  the  relativistic  core  orbitals  of  Pb 
arc  less  diffuse  than  the  analogous  nonrelativistic  ones,  the  relativistic  4f,4f*.  5d  and  Sd*  are  slightly  more  dif¬ 
fuse.  Thercfote  the  relativistic  contraction  of  the  5s,  and  5p  is  largely  offset  by  the  4f  and  5d  functions,  particu¬ 
larly  in  the  local  potential.  Qualitative  predictions,  such  as  the  occurrence  of  shorter  bond  lengths  in  the  relativistic 
case,  based  on  the  assumption  that  all  relativistic  core  orbitals  are  more  contracted,  arc  not  necessarily  valid. 

It  is  also  useful  to  compare  these  results  with  those  of  Schwenzer  et  al.  [2),  who  performed  a  unique  minimal- 
basis  set  all-electron  SCF  calculation  on  PbO.  These  authors  compared  their  results  with  experiment  in  order  to 
estimate  the  magnitude  of  relativistic  effects.  It  is  interesting  to  note  that  our  nonrelativistic  results  are  in  appre¬ 
ciably  better  agreement  with  experiment  than  theirs  (1.871  A  for  the  bond  length  and  868  cm  1  for  the  frequen¬ 
cy),  indicating  that  most  of  their  discrepancy  with  experiment  was  due  to  the  limited  basis  set  employed,  rather 
than  to  use  of  the  nonrelativistic  hamiltonian.  , 

We  conclude  that  for  the  observable  properties  in  table  1  the  relativistic  and  nonrelativistic  results  are  nearly 
the  same.  This  is  in  agreement  with  the  general  conclusions  of  Schwenzer  et  al.  but  not  with  those  of  ,Desclaux  and 

* 

*  The  oxygen  orbitals  are  from  ref.  (8).  The  lead  orbitals  and  basis  were  obtilned  from  atom  optimisation.  Details  of  the  lead  basis 
set  will  be  published  elsewhere.  *- 

\ 

Table  I  \ 

Computed  properties  of  PbO  _  _ 


Nonrelativistic 

Relativistic 

Experimental 

orbital  energies  (au) 

*|)0 

-1.120 

-1.274 

*140 

-0.464 

-0.878 

- 

*j  70 

-0.283 

-0.418 

- 

*9* 

-0.396 

-0.484 

- 

bond  length  (A) 

1.910 

1.903 

1.922 

vibrational  frequency  (cm'1 

)  749 

760 

721 
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Pyykko  1 1 1,  wlio  found  a  much  larger  difference  in  Pbll^  using  the  spherical  approximation  The  results  found  by 
Das  and  Wahl  1 10)  using  a  different  valence-electron  method  suggest  that  relativistic  effects  may  be  much  more 
import. ml  in  properties  such  as  electronic  transition  cnetgtes,  rather  than  in  ground-state  properties. 
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Abst  ract 

A  general  theoretical  approach  and  computational  method  are  described 
tor  employing  pseudopotential  operators  to  simplify  ab  initio  molecular  multi¬ 
configuration  SCF  calculations.  Considerable  computational  effort  may  be  saved 
in  this  way  since  only  the  valence  orbitals  need  be  computed  explicitly.  Using 
the  ground  electronic  state  of  l.i  ,  as  an  example,  the  role  of  the  I.agrange  mul¬ 
tipliers  that  arise  in  the  MCSCF  formalism  is  discussed.  The  pseudopotential 
and  effective  potential  approaches  described  herein  are  shown  to  accurately  reproduce 
the  detailed  electronic  structure,  including  configuration-mixing  coefficients, 
of  Li,. 


i 


1.  Introduction 


In  recent  years  there  has  been  an  increasing  interest  in  techniques  for 
carrying  out  ab  initio  studies  of  atomic  and  molecular  properties  treating 
only  valence  orbitals  explicitly.  Such  techniques  are  greatly  extend  the 
number  of  chemical  systems  that  are  amenable  to  precise  theoretical  calculations. 
In  particular  they  facilitate  the  study  of  molecules  containing  heavy 
elements,  in  which  the  numbers  of  core  electrons  is  large.  But  in  principle  the 
valence-only  approach  allows  most  of  the  properties  of  any  systenj,  for  a  given 
amount  of  computational  effort,  to  be  found  more  precisely. 

We  have  previously  described  a  general  theoretical  approach  that  allows 

carrying  out  sel f-consistent- field  calculations,  treating  only  valence  orbitals, 

1  2 

that  parallels  conventional  SCF  methods,  ’  We  now  wish  to  extend  our  theo¬ 
retical  development  to  the  more  general  cases  of  open-shell  SCF  and  multicon¬ 
figuration  SCF  (MCSCF)  applications.  This  extension  is  clearly  necessary  in 

order  to  predict  important  properties  -  for  example,  dissociation  energies  of 
inorganic  compounds. 

Other  workers  have  reported  successful  valence  molecular-orbital  cal¬ 
culations  using  open-shell  SCF  and  MCSCF  methods . ^  6  These  have  all  been 
based  on  what  might  be  termed  the  core-potential  approximation,  i.e.'  that  the 
effect  of  the  core  electron  distribution  on  the  valence  may  be  represented  by 
potential  functions  characteristic  of  the  core  but  independent  of  the  valence 
orbitals.  This  approximation  simplifies  calculations,  but  its  accuracy  and 
theoretical  basis  have  not  been  studied.  1 

In  this  paper  we  will  show  how  one-electron  Phi 1 1 ips-Kleinman  pseudo- 
potentia!  operators,  which  replace  core-valence  orthogonality  constraints  in 
valence-only  calculations,  may  be  employed  in  MCSCF  calculations. 

This  derivation  is  nonempirical  and  follows  from  the  variation  condition  im¬ 
posed  on  the  total  energy.  The  sum  of  all  terms  that  simulate  the  effect  of 


core  orbitals  on  the  valence  we  refer  to  as  an  effective  potential  to  distin¬ 
guish  it  from  pseudopotentials  (although  the  former  contains  nonlocal  operators 
and  hence  is  not  a  potential  function  in  the  usual  sense).  The  terms  in  addi¬ 
tion  to  the  pseudopotential,  representing  electrostatic  interactions,  have  been 
derived  in  an  earlier  paper.  “  We  will  therefore  briefly  review  their  proper¬ 
ties  and  show  how  their  definition  follows  from  the  form  of  the  pseudopotential. 

The  computational  method  based  on  this  theoretical  approach  allows 
parol  lei  all-electron  and  valence-only  calculations  to  be  carried  out  and  com¬ 
pared.  We  report  such  a  study  using  one-configuration  descriptions  of  Li  and 

I. i,  and  a  five-configuration  description  of  L^-  These  results  demonstrate  both 
the  validity  of  the  individual  approximations  involved  in  using  pseudopotential 
operators,  and  also  the  precision  of  the  overall  method. 

II.  General  Theory 

The  greatest  simplification  in  the  description  of  the  valence-electron 
distribution  results  if  the  valence  orbitals  are  allowed  to  be  completely  inde¬ 
pendent  of  the  core,  and  hence  are  not  required  to  have  nodes  in  the  core  re¬ 
gion.  Therefore  in  this  section  we  will  discuss  the  conditions  that  must  be 
met  by  a  set  of  nodeless  spacial  functions  (pseudoorbitals)  in  oruer  that  they 
be  numerically  equal  to  the  equivalent  valence  canonical  molecular  orbitals  in 
the  valence  region.  Here  valence  will  be  defined  as  all  points  in  space  where 
the  sum  of  all  core  orbitals  is  negligibly  small.  Note  that  these  conditions 
alone,  rather  than  consideration  of  computed  properties,  lead  to  the  formation 
of  a  pseudopotential.  For  example,  a  parameterization  scheme  based  on  atomic 
orbital  energies  might  precisely  reproduce  molecular-orbital  energies.  But 
all  the  properties  of  a  system  may  in  principle  be  obtained  from  its  wavetunc- 
tion;  and  such  parameters  giving  accurate  orbital  energies  may  well  not  suffice 
for  a  good  estimation  of  some  physically  important  property.  Also,  it  is  nearly 
always  found  that  when  the  pseudoorbitals  are  properly 


28 


determined  they  are  sufficiently  similar  to  the  true  orbitals  that  many  other 
properties,  such  as  molecular  geometries,  are  also  accurately  reproduced. 

The  condition  that  the  total  energy  obtained  by  any  SCF  method  be  a 

•  i  7 

minimum  with  respect  to  spacial  variations  of  an  orbital  ^  may  be  written 

*  I  <!)>  .  |  =0  (1) 

jA  J  k  J 

where  l-k  is  the  Fock  operator  for  orbital  k.  The  terms  <4> .  |  Fk  1 4>,  >  are  Lagrange 

J  k 

multipliers  that  ensure  orbital  normalization  when  j  =  k  and  orthogonality  for 
J  ^  k .  In  these  expressions  and  those  below  we  will  use  indices  j  and  k  to  re¬ 
fer  to  orbitals  ordered  regardless  of  symmetry,  and  m  and  n  within  a  symmetry 
block . 

In  the  MCSCF  method,  the  wavefunction  is  also  expanded  as  a  superposi- 
„  8,9 

tion  of  electronic  cont igurat ions .  These  are  generally  obtained  by 

one  or  a  pair  of  electrons  from  the  Hartree-Fock  configuration  to  excited  orbi¬ 
tals,  each  of  which  is  optimal  in  the  sense  that  it  satisfies  its  own  Fock  equa¬ 
tion  as  in  Eq.  (1).  Excitations  are  chosen  as  connecting  valence  orbitals, 
since  electron  correlation  in  the  region  near  the  nuclei  is  of  little  practical 
interest.  In  the  total  energy  found  by  this  procedure,  all  terms  depending 
solely  on  the  core  orbitals  appear  as  ar.  additive  constant  as  the  valence  orbi¬ 
tals  are  varied.  There  is  thus  a  logical  relationship  between  the  molecular 
MCSCF  approach  and  its  valence-only  counterpart. 

The  condition  that  must  be  satisfied  in  order  that  the  energy  be  a 

minimum  with  respect  to  variation  of  the  contributions  of  all  configurations 

9 

in  the  wavefunction  may  be  written  in  the  form 

7  A  H  =  EA 

*■  no  mo.no  mo 

n  9 


(21 
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whore  I;  is  the  total  enciey  and  A  is  ’he  vector  of  configuration  expansion 
coefficients  of  the  mth  root  of  li,  assumed  to  be  the  lowest,  of  symmetry  o.  In 
forming  !1  one  may  consider  either  the  requisite  integrals  over  all  the  orbitals, 
or  equivalently  those  which  involve  only  valence  orbitals.  The  values  of 
may  be  shown  to  be  the  same  in  either  case;  but,  in  the  latter,  li  -  1:  the  total 

energy  minus  the  self-energy  of  the  core.  Assume  for  simplicity  that  each  con¬ 
figuration  may  be  associated  with  either  the  11a rt ree- Fock  case  or  an  excitation 
to  a  specific  higher  orbital.  Then 
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o  m ,  o  m 


N  (h 
om  am 


In, 


ii* 


p  ) 

pi  am, pi 


d  K 

o  am, am 


( 3a ) 


11  =  (.d  d  )  K 

ump  n  o  p  om.pn 


for  am  /  pn, 


f  3b ) 


where  N  is  the  occupation  number  of  orbital  m,  h  is  the  one-electron  Hamil- 
am  r  am 

tonian,  dQ  is  the  dimension  of  symmetry  o  (i.e.  1  for  sigma,  2  for  pi,  etc.)  and 


Am, pi 


(d.d  ) 
A  p 


U 


Am, pH 


^Am.pJ 


(4) 


Here  J  and  K  .  are  the  conventional  Coulomb  and  exchange  terms  connecting 

Am, pi  Am, pi 

the  valence  orbitals  with  each  other  and  with  the  core.  The  solutions  of  (1)  and 
(2)  are  interdependent,  since  the  matrix  elements  of  H  depend  on  the  spacial  form 
of  the  orbitals,  and  each  Fk  may  be  shown  to  depend  on  the  A^ . 

In  constructing  the  H  matrix  in  a  valence-only  calculation,  it  should 
be  noted  that  any  effective  potential  employed  to  simulate  the  core  or¬ 
bitals  is  required  to  be  a  one-electron  operator,  so  that  it  modifies  only 

H  in  (3).  The  off-diagonal  terms  H  depend  solely  on  the  (two-electron) 

om,om  om,pn 

exchange  integral  K  .  Therefore  A.  ,  the  expansion  of  the  wavefunction  in 
6  °  om,pn  10 
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terms  of  configurations,  is  independent  of  the  effective  potential,  except  for 

the  indirect  dependence  through  Lq .  (1).  Any  calculation  employing  one-electron 

operators  to  represent  the  core  orbitals  must  rely  on  the  assumptions  that  (i) 

the  solutions  of  to.  (1)  are  not  sensitive  to  the  values  of  A.  ,  and  (ii)  the 

10 

pseudoorbi ta 1 s  obtained  in  this  way  are  equal  in  value  to  the  true  orbitals  over 
a  sufficiently  large  region  of  space,  so  that  two-electron  integrals  over  both 
sets  of  functions  are  approximately  the  same. 

The  form  of  the  pseudopotential  and  effective  potential  must  therefore 
be  found  by  relying  on  Eq.  (1).  Let  xk  he  the  pseudoorbital  derived  from  4>k  by 
adding  an  arbitrary  linear  combination  of  core  orbitals.  If  P  is  the  core-pro¬ 
lection  operator. 


P 


I 

j  =core 


1 j > <4>  j  I 


(5) 


then 


♦k  =  (1  *  P)*k- 


(bj 


Dividing  the  <f>m  into  core  and  valence  functions,  Eq.  (1)  may  be  written  in  the 
form 


F'C^lt  "  1  <(<>  1  I  I  1  "  £  I 

j=core  J  J  H=valenceA 

=  <<*>k I  pK I <*’k>4>k *  (7) 


Partially  substituting  (6)  into  (7), 
„k 


Fk(l  -  PHk  -  l  «UFk|(l  -  P)xk>*:  -  l  «frn|FklH>*t 
j=core  J  J  t=valerice^k 


=  <*JFk  |  ♦•>  (1  -  P)xt- 


(«) 
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Now  it'  wo  lot  represent  ^  4.  L  1 1;  ;  4>  k  >  ,  the  diagonal  Lagrange  multiplier,  and  use 
definition  (  ;>)  , 


rkll  -  l’)xk  -  pi  k  1 1 1  -  P)xkl  *  l  '♦i|Fkln>H  *  nkc  1  -  l*ixk 

J-valenco^K 


C.»  ) 


After  rearranging,  tliis  becomes 


t‘kx,  ♦  ( PFk I*  -  f'kP  -  l’l-K  ♦  ftk  1 x k  "  £  ^  i 


,.k„  ,.k 


•k, 


£- va  1  ence/fk 


*kx  k 


in  I 


Hie  term  in  brackets  is  now  the  usual  form  of  the  generalized  Phi  1 1 ips- kleinman 

pseudopotential  for  orbital  k.  It  may  be  reduced  to  an  expression  involving 

only  overlap  integrals  if  it  is  assumed  that  the  core  orbitals  are  approximate 

k  k 

eigenfunctions  of  each  F  ,  so  that  P  and  F  commute.  In  this  case  the  pseudo¬ 
potential  operator  becomes 

,PP,  ,  .  „ 


V  (n.)  i  n  P  -  FkP  =  l  K  -  nJU  •><*■! 

j=core  )  )  ) 


I  1  1  i 


This  now  allows  us  to  form  an  equation  for  analogous  to  bq.  (1).  If  we  define 


•gk  _  r-k  l/PP/  ■> 
F  =  F  +  V  (nk J 


(12) 


then  Eq.  (10)  may  be  written  in  a  form  in  which  core-valence  Lagrange  multipliers 
do  not  appear, 


Xk  *  l  <xjFk|xk>X, 

i=valence^k 


Vk 


(  1  s  I 


providing 


l  <^0 1  pK  I <t'k>4>fi.  £  ^  <x«. !  pkj xk>x 

£  =  valerice^k 


Q  1  *  I  A  1  A  Q 

t=valence^k 
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'♦klF  !♦. 


!  i  i 

'k1  1  'k  ’ 


(14b) 


K'vse  conditions  arc  genera  ]  1  y  satisfied  in  practice  within  a  good  approximation. 
1-or  example,  all  the  terms  appearing  i  n  (11a)  are  generally  quite  small,  and 
vanish  in  the  case-  ot  one- con t l gural ion  methods.  Also  since  valence-only  calcu¬ 
lations  generally  reproduce  valence  orbital  energies  quantitatively  in  one- 
conf l gurat ion  calculations,  (I4hj  is  expected  to  hold  at  least  tor  orbitals  that 
are  not  highly  excited. 

In  applying  these  equations,  and  n .  for  core  orbitals  may  be  taken 
from  readily  available  reference  systems,  such  as  the  constituent  atoms  of  a 
molecule.  The  and  are  then  found  hy  solving  fcq .  (13) . 

Note  that  each  orbital,  including  those  not  occupied  in  the  Hartree- 
lock  configuration,  has  its  own  pseudopotential  operator.  In  conventional  MCSfF 
calculations  using  a  small  number  of  configurations,  all  orbitals  that  are  com¬ 
pletely  filled  in  every  configuration  share  a  common  Fock  operator.  Pseudopo- 
tential  calculations,  on  the  other  hand,  require  a  distinct  F  operator  for 
each  occupied  orbital,  regardless  of  its  degree  of  occupation.  This  does  not 
pose  a  computational  difficulty,  however,  and  the  difference  in  the  number  of 
these  operators  disappears  when  the  number  of  configurations  is  large. 

Eq.  (13)  is  much  more  readily  solved  than  tq.  (1)  since  fewer  functions 

need  to  be  found  and  these  may  be  described  by  a  much  smaller  set  of  basis  func- 

tions.  However,  F  contains  a  large  number  of  two-electron  integrals  connecting 

core  and  valence  orbitals.  In  order  to  solve  (13)  it  is  useful  to  introduce  a 

set  of  potential  functions,  expressible  m  terms  of  one-electron  integrals,  to 

replace  them.  We  have  previously  described  a  technique  for  generating  these 

2 

potentials  for  one-configurat ion  calculations.  “  The  same  formulas  apply  in 
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the  inul t iconf igurut ion  case  since,  assuming  all  core  orbitals  are  closed,  pre¬ 
cisely  the  same  type  of  core-valence  Coulomb  and  exchange  integrals  occur. ^ 

liq.  (13)  must  therefore  be  reformulated  in  terms  of  a  set  of  Fock  op¬ 
erators  that  depend  on  the  coordinates  of  only  valence  electrons.  In  general 
F  is  to  be  replaced  by  +  V* ^(n^)  +  ,  where  is  a  one-electron  operator 

whose  form  is  to  be  determined.  It  is  again  necessary  to  choose  a  reference 

L 

system,  generally  the  free  atoms,  and  to  assume  that  once  W  is  found  for  that 

case  it  will  be  t ransferrable  to  other  systems,  such  as  molecules. 

Note  that  in  Eq.  (13)  the  summation  vanishes  by  symmetry  for  a  one-conf i gurat ioi 
description  of  an  atom. 

lit  ere  are  thus  three  possible  ways  to  define  W  :  require  that  it  (i)  reproduce 

.  — k  .  .  4 

the  form  of  F  ,  (ii)  reproduce  the  value  of  n^,  or  (iii)  minimize  (F  - 

Ihe  second  choice  may  be  the  most  accurate  in  some  cases,  since  it  can  compen- 

k 

sate  for  lack  of  commutivity  between  F  and  the  core  orbitals,  and  it  causes 
condition  (14b)  to  be  satisfied  by  construction.  However  this  method  relies  on 
only  a  small  set  of  distinct  numbers,  atomic-orbital  energies,  and  hence  does 
not  allow  a  very  detailed  description  of  W  .  Choosing  its  functional  form  care¬ 
fully  gives  good  results  in  many  cases,  but  in  general  only  when  small  basis 

I 

sets  are  employed?  The  third  choice  above  assumes  that  Eq .  (13)  can  be  satis¬ 
fied  at  all  points  in  space.  This  Ciinnot  be  true  if  x^  is  described  by  a  finite 
linear  combination  of  basis  functions.  In  other  words,  such  a  definition  makes 

a  sensitive  function  of  x^’  Further,  if  greatly  differing  descriptions  of 

k 

x ^  are  used  in  the  atom  and  molecule,  W  so  determined  for  an  atom  may  not  be 
appropriate  to  the  molecular  case. 

Therefore,  as  in  our  earlier  work,  we  impose  the  operator-equivalence 

-  k  pp  k  — k 

|F  ♦  V  (n^)  ♦  W  |xk  =  F  xk  (lb) 


;ondi t ion 


which,  from  f-.q .  (10),  may  he  written 
-  k  k 

U'  +  W  )xk  -  l'xk.  (lb) 

This  is  clearly  the  same  as  minimising  ( I-*'  -  'ik)\k  for  1  a  rye  valence  basic  sets. 
It  will  also  reproduce  nk  provided  Eq.  I  14b)  is  satisfied,  although  of  course 
the  converse  is  not  true. 

An  additional  advantage  of  defining  Wk  in  operator  form  is  that  Eq . 

(  H>  I  can  generate  a  potential  commensurate  with  any  function,  rather  than  only 
with  a  specific  valence  atomic  orbital.  Suppose  for  example  that  xk  is  described 
by  a  series  of  basis  functions,  ij/.  , 

Xk  =  I  c\  (17) 

1 

There  are  two  ways  of  determining  W,k ,  which  are  equivalent  for  the  ground  state 
of  an  atom  but  not  otherwise.  First,  W*"  may  be  associated  with  a  particular 
atomic  angular-momentum  type,  such  as  s  or  p.  Then  Eq .  lib)  may  be  solved  set¬ 
ting  xk  equal  to  the  atomic  pseudoorbital  of  each  type.  More  generally,  W*^ 
may  also  be  determined  separately  for  each  member  of  the  basis  set, 

W-  »  UFk  -  Fk)*i|/*i.  (IK) 

Using  (17)  this  gives 


In  other  words  W  may  be  considered  as  being  projected  either  onto  each  angular- 
momentum  type  or  onto  each  basis  function,  both  alternatives  satisfying  l.q .  (h<) 


for  atoms. 
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F^  and  are  of  course  independent  of  k  in  the  closed  shell  Hartree- 
Fock  case.  In  general  all  the  1^  (and  F*' )  operators  are  quite  similar  except 
tor  highly  excited  valence  orbitals.  The  difference  ( F^  -  F^)  is  essentially 
independent  of  k  in  all  cases,  since  it  reflects  only  core  electron  density, 
he  will  therefore  drop  the  k  in  l.q .  (IS). 

lor  any  other  valence  function  x^.  which  may  or  may  not  approximate  an 
atomic  orbital,  we  have 


IF*  ♦  W*  ) x , 


:X 


F'x .  ♦  l  h  C  V  =  v\o  *  W  -  lX)C.\. 

t  Fill  *8.  F  11 

i  l 


F\.  (-'<•)) 


Thus  the  core  potential  determined  in  this  way  would  be  expected  to  reproduce  to 
some  extent  the  behavior  of  excited  molecular  orbitals  such  as  those  formed  from 
excited  or  ionized  states  of  the  constituent  atom,  even  though  is  generally 
based  solely  on  ground  states  of  neutral  atoms. 

Other  features  of  these  core  potential  functions  needed  to  carry  out 
molecular  calculations  are  the  same  as  have  been  described  previously. 

Finally,  note  that  specifying  the  form  of  the  lock  operator  is  equiva¬ 
lent  to  setting  forth  an  expression  for  the  total  energy.  Therefore  for  consis¬ 
tency  the  effective  potential  should  modify  the  other  variation  condition  of  the  energy, 

pp 

Fq.  (.-),  as  well.  It  may  he  readily  seen  that  this  simply  involves  adding  V  In  i  to 
h  and  replacing  P  „  bv  matrix  elements  of  W.  for  all  core  values  of  at  in  Fq .  (3a), 

III.  Representative  Calculations  on  Li, 

In  this  section  we  will  examine  the  accuracy  of  some  of  the  approxima¬ 
tions  made  above,  and  illustrate  the  type  of  results  that  are  obtained.  These 
will  he  based  on  the  ground  electronic  state  of  Li-,,  the  simplest  molecule 
possessing  core  orbitals.  It  was  also  one  of  the  first  molecules  to  he  exten¬ 


sively  studied,  by  lias,  using  the  NK'.SCI  method.’  Throughout  we  will  assume  an 
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°  1  ] 

o xj>e r i inont ;i  1  value  u f  2.(i73  A  for  the  equilibrium  '  >nd  length.  Our  basis  set 

is  the  molecularly-optimized  2-t;  set  recently  developed  by  Dixon,  Gale  and 

Jordan .  lor  consistency  the  s- symmetry  components  of  this  basis  set  are  used 

) 

to  describe  the  reference  system,  the  ground  “S  state  of  the  free  Li  atom. 

The  W.  functions  were  generated  for  each  valence  basis  function  as 
described  previously,  using  liq.  (JO)  of  Ref.  2.  The  single  atomic  pseudoorbi¬ 
tal  is  just  the  valence  part  of  the  canonical  atomic  orbital,  that  is  only  the 
part  that  depends  on  2s  basis  functions.  The  core  component  of  each  is 
fixed  by  this  definition.  Note'  that  by  using  the  method  described  above  it  is 
not  necessary  to  postulate  the  form  of  a  Li  2p  orbital  in  order  to  form  any  W. . 
for  carrying  out  computations  it  is  convenient  to  expand  each  function  in  the 
fo  rm 


W  •  (r) 
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where  r  is  the  distance  from  the  nucleus.  The  four  resulting  potential  func¬ 


tions  in  terms  of  d  ,  and  a .  .  for  each  function  ip  -  are  shown  in  Table  I. 

l  i  li  l 

Our  multiconfiguration  molecular  results  are  based  on  the  five  (six 


real)  configurations  that  are  most  significant  in  terms  of  both  participation 


in  the  wavefunction  and  minimization  of  the  total  energy.  These  are 
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where  clearly  la  and  lo  represent  the  core  orbitals.  2a  and  2o  are  formed 
8  “  g  u 

primarily  from  Li  2s  atomic  orbitals,  In  from  2p,  and  3o  and  4a  from  2p  and 

u  8  g 

3s.  At  infinite  internuclear  sepa rat i on,  t h is  set  of  configurations  reduces  to 

TO")  "7  T  7 

the  Ha  rt  ree- Lock  description  of  the  free  atoms,  i.e.  1 //2  ( lo“  lo  ^  2o“  -  lo~lo  do''). 

8  u  g  g  u  u 

In  forming  the  II  matrix  we  have  included  all  terms,  including  intra- 

atomic  core-core  interactions,  in  order  to  allow  comparison  of  L  with  H  as  descri 

a 

The  basic  assumption  in  deriving  the  equation  for  the  effective  poten¬ 
tial  has  been  that  the  core  is  not  strongly  affected  by  differences  in  its  en¬ 
vironment  ot  valence  orbitals.  Specifically  we  require  that,  to  a  good 
approximation,  the  core  orbitals  are  the  same  in  the  atom  as  in  the  molecule  so 

_  L 

that  they  may  be  used  to  construct  accurately  each  F  operator.  Further  we 

must  postulate  the  existence  of  one  set  of  core  functions  that  commutes  with 
.k 

each  of  the  F  .  (In  Li.,  the  last  condition  does  not  apply  to  it  orbitals.)  A 
useful  measure  of  the  extent  to  which  all  these  requirements  are  met  is  compar¬ 
ison  of  the  diagonal  Lagrange  multipliers,  n^.  •  As  a  test  of  the  first  condition, 
we  list  in  Table  II  the  values  of  for  the  core  orbitals  obtained  by  a  conven¬ 
tional  SCF  calculation  on  Li,  and  the  commensurate  single-configuration  and 
multiconfiguration  results  for  Li^.  As  might  be  expected,  all  these  values  are 
about  the  same  -  the  major  differences  being  between  the  atomic  and  the  two 
molecular  calculations. 

The  general  similarity  of  the  core  solutions  obtainable  from  valence 
F^  operators,  as  required  by  our  second  condition  above,  is  well  known.1'  But 


we  must  also  require  that  atomic  core  orbitals  are  approximate  eigenfunctions 

g 

of  each  F  .  The  latter  might  be  referred  to  as  virtual  core  orbitals,  since 
within  the  SCF  form. il ism  they  are  not  allowed  to  be  occupied.  The  eigenvalues 
corresponding  to  these  core  lo  -  and  lo^-type  eigenfunctions  for  the  Fuck 
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operators  of  each  of  the  valence  orbitals  in  the  mul  t  iconf  igurat  ion  calculation 

are  shown  in  Table  111.  (The  11  lock  operators  possess  no  such  solutions.)  The 

results  tor  the  Jo  orbital  are  very  close  to  those  obtained  from  the  lo  and 
g  g 

lo  shown  m  Table  11.  These  three  would  be  identical  in  a  one-configuration 
calculation.  They  are  again  quite  close  to  the  atomic  value.  (he  higher  orbi¬ 
tals  show  a  progressive  decrease  in  their  core  eigenvalues  as  they  become  less 

like  atomic  Js  orbitals.  This  difference  is  small  except  for  4o  .  As  shown 

g 

below,  the  4o^  orbital  contributes  relatively  little  to  the  molecular  wave- 
function,  so  that  the  ability  of  a  pseudopotential  to  describe  it  as  accurately 
as  the  others  is  not  critical.  'Hie  same  would  also  be  true  of  a  larger  config¬ 
uration  list,  which  would  include  primarily  other  occupations  of  the  first  four 

valence  orbitals  (2o  ,  2«j  ,  In  ,  and  3o  ). 

g  u’  u’  g' 

I’P 

The  one-electron  form  of  the  pseudopotential  operator,  V  (n^)  as  de¬ 
fined  by  Eq .  (11)  depends  on  the  differences  between  core  and  valence  Lagrange 
multipliers.  In  Table  IV  we  list  the  values  of  for  the  valence  orbitals  of 
Li,  and  one-  and  five-configuration  descriptions  of  Li2-  Die  atomic  Js  and 
molecular  2a ^  values  are  about  the  same,  but  all  those  corresponding  to  excited 
orbitals  are  again  found  to  be  more  negative.  Comparing  values  of  Table  III  with 
those  from  the  last  column  of  Table  IV,  there  is  found  to  be  a  roughly  parallel  de¬ 
crease  in  n  for  both  core  and  valence  u  orbitals.  This  suggests  that  an  approximate 
pseudopotential 

based  on  atomic  core  Lagrange  multipliers  but  setting  =  n 7  of  Li  for  all  k 
would  be  reasonably  accurate  in  this  case,  since  it  would  roughly  reproduce  the 
differences  in  seen  in  the  all-electron  calculation. 

We  can  also  use  the  values  of  for  valence-orbitals  to  check  condi¬ 
tion  (14b)  in  conjunction  with  the  W.  functions  by  carrying  out  parallel  all- 
electron  and  valence-electron  calculations.  In  the  latter  case  we  employ  Eq . 

(13), add  the  analogous  effective  potential  to  h  in  Eq .  (3a),  and  define  the 


configuration  list  to  be  just  2o^,  2ou>  etc.  The  results  are  shown  in  Table  V, 
along  with  the  core  projection  of  each  valence  pseudoorbital,  defined  as 


I  <xk 
j  =core 


which  measures  the  extent  of  its  interaction  with  the  core.  The  agreement  with 
the  results  in  Table  IV  is  good  lor  the  2o^  orbital  but  is  again  somewhat  less 
satisfactory  for  the  excited  orbitals.  For  3o  and  4o  ,  however,  the  overlap 
with  the  core  is  relatively  small  so  that  F^  is  less  sensitive  to  for  these 
orbitals . 

We  now  wish  to  examine  the  solution  of  Eq.  (2)  in  terms  of  this  valence- 
electron  description.  In  particular  we  need  to  varify  that  the  exchange  inte¬ 
grals  in  Eq.  (3)  calculated  over  pseudoorbitals  form  a  valid  approximation  to 
those  obtained  using  the  true  orbitals.  Table  VI  shows  the  configuration-inter¬ 
action  matrix  H  obtained  in  these  two  cases.  Clearly  the  off-diagonal  elements 
defined  by  Eq.  (3)  are  reproduced  with  reasonable  accuracy.  The  two  matrices 
are  seen  to  differ  by  essentially  a  constant  times  the  identity  matrix.  This 
constant  represents  the  self-energy  of  the  core  electrons,  E  -  E^  in  the  all¬ 
electron  case.  It  follows  that  the  eigenvectors  describing  the  contribution 
of  each  configuration  to  the  total  wavefunction  should  also  be  approximately 

the  same.  This  is  demonstrated  in  Table  VII.  The  value  of  A  coefficients 

mo 

obtained  from  the  valence-electron  calculation  are  quite  accurate  for  each  or¬ 
bital,  including  4a^,  which  however  makes  only  a  small  contribution  (.007)  to 
the  wavefunction. 

Finally,  it  is  useful  to  compare  the  total  atomic  and  molecular  ener¬ 
gies  obtained  using  an  effective  potential  with  those  obtained  from  the  equi¬ 
valent  all-electron  method.  These  values  as  shown  in  Table  VIII  for  the  Li 
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atom  aiul  the  one  -  con f i gurat ion  and  five-configuration  Li,  calculations.  Compar¬ 
ing  the  second  and  third  columns,  the  valence-only  results  somewhat  overestimates 
the  correlation  energy,  giving  .OS-17  a.u.  compared  to  .02t>8  a.u.  in  the 
all-electron  case.  However  comparing  the  first  and  third  columns  shows  that  the 
dissociation  energy  is  quite  well  reproduced,  being  .0246  in  the  former  case  and 
.U2t)4  in  the  latter.  The  experimental  value  is  .0377.  Due  to  the  small  number 
of  configurations  we  have  employed,  our  result  is  also  less  than  that  obtained 
by  the  large  configurut ion- interact  ion  calculation  in  Kef.  12.  These  authors, 
using  the  same  basis  set,  obtain  a  dissociation  energy  of  .0320  a.u. 

IV.  Conclusions 

In  the  two  preceding  sections  we  have  (a)  outlined  a  theoretical  approach 

for  computing  valence  molecular  orbitals,  (.bj  demonstrated  that  the  approximations 

involved  are  valid,  and  (c)  shown  the  type  of  results  that  may  he  obtained  by  its 

use.  Since  at  present  molecular  calculations  are  most  readily  carried  out  by 

the  expansion  method,  it  is  useful  that  our  approach  may  be  carried  out  entirely 

within  a  particular  basis  set.  This  obviates  the  artificial  adjustments  that 

arise  in  attempting  to  interface  an  analytical  molecular  wavefunction  with  a 

numerical  (Hartroe-Fock  limit)  atomic  calculation,  as  in  the  method  of  Kahn  and 
14 

Goddard.  Note  also  that  the  calculation  of  dissociation  energies  using  the 
method  described  above  does  not  require  consideration  of  a  potential  function 
"tail",***  which  is  an  artifact  of  their  method.  Since  we  have  not  vet  optimized 
or  parameterized  our  procedure  in  any  way,  it  seems  likely  that  doing  so  could 
improve  its  accuracy  further.  However  we  feel  that  such  parameterizing  is  gen¬ 
erally  unreliable,  and  would  not  give  constructive  results  in  the  long  run. 
Finally,  it  should  be  emphasized  that  our  overall  valence-orbital  procedure  is 
sufficiently  efficient  that  computations  on  large  molecular  systems  are  feasible 
using  this  method  in  its  present  form. 


Table  I.  Seven-term  expansions  ot'  the  core-valence  interaction  potentials 
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Table  11.  Core  diagonal  Lagrange  multipliers  of  Li  (Is) 
and  Li^  (lo  ,  lo^)  in  the  SCF  and  MCSCF  approximations 
(a.u. ) . 

Li 

1  configuration  -2.4772 

l«2 

1  configuration  -2.4449,  -2.4446 
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Table  111.  Virtual  core  eigenvalues  of  the  four  a - 
symmetry  Lock  operators  in  a  fi  ve-  confi  gurat  ion  descrip- 
t ion  of  Li  ,  (a . u. ) . 


I’ock  operator 


2o 

u 


3o 

U 


•to 

8 


core  eigenvalues 


lo 

u 


-2.4459,  -2.4455 
-2.4783,  -2.4781 
-2.5142,  -2.5135 


-3.2480,  -3.1778 
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Table  V.  Valence  Lagrange  multipliers,  n^,  for  each 
valence  pseudoorbital  found  from  an  effective -potential 
calculation  on  l.i  ,  using  five  optimized  valence  config¬ 
urations.  Also  shown  is  the  projection  of  each  function 
on  the  core  orbitals,  as  defined  in  the  text. 


orbital  core  projection 


2o  0.0447 

& 

2o  0.0644 

u 

3o  0.0141 

6 

4a  0.0293 

8 

In 

u 


-0.1934 
-0.2980 
-0 .2987 
-0.4825 


0 


-0.2577 
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Table  VI.  l.'nnf  igurat  inn  -  interact  ion  matrices  resulting  from  a  five- 
configuration  desc  ri  1 1 1  1  on  of  Li  ,  in  the  valence  electron  mu.  a  1 1 -elect  roll 
approx  niiat  ions .  The  columns  of  eacli  matrix  are  ordered  in  the  same  way  as 
the  rows.  Since  these  matrices  are  symmetric,  only  the  bottom  half  is 
shown . 


conf i gur at  ion  valence-electron  case 


» 

jj  ~  -o..viyr> 


1  r 

U 

0. 0-108 

-0.2285 

> 

Jj  " 
u 

0 . 0440 

0.021  7 

-0.1524 

i 

3o  " 

K 

i 

4o 

8 

0.0455 

0.0151 

0.0087 

0.0085 

0.0188 

0.0313 

-0.1582 

0.0156 

0.1787 

conf i gurat ion 

al  1 

-electron  case 

[core)  2a  ^ 

-14.8288 

■) 

[core]  l"u~ 

0.0414 

-14.6504 

[core]  2o  ^ 

0.0448 

0.0220 

-14.5809 

[core]  3o ^ 

0 . 04  30 

0.008b 

0.0160 

-14.6082 

[core]  4o^ 

0.0143 

0.0083 

0.0314 

0.0133 

-14.2847 

t 


i 
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1 .»!» 1 Ul.  Cmili  yur.it  i  on  expansion  coefficients  of  the 


(■rmihii  c  1  cc  t  run  i  c 

e  1  ci  t  roll  aiul  a  1  1  -e 

state  of  Li  ,  resnl t l ny 

lec tron  calculations. 

from  valence 

con t  i^i i rat  ion 

va  1  ence-e  1  ect  roll 

al  1  -  electron 

(cure|  Jo 

0.9533 

0.9S2S 

| Corel  In 

11 

> 

-0. 2075 

-0 . 169S 

I corel  Jo 

u 

> 

-0.1279 

-0.1187 

1 corel  So 

& 

-0. 162b 

-0.1454 

[corel  4  a 

& 

-0.0072 

-0.0091 
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Relativistic  Core  Binding  Energies  of  Selected  Atoms:  Comparison  with 
Experiment  and  Other  Calculations 

R.  J.  Key,  M.  S.  Banna  and  C.  S.  Ewig 

Single-configuration  Dirac- lock  binding  energies  are  reported  for  K  2p 
kb  3p ,  Cs  3d,  Mg  Is,  Zn  2p,  and  Cd  3d  a-  well  as  for  the  corresponding  levels  i 
cations  of  these  elements  and  the  neighboring  rare  gases.  The  results  are  com¬ 
pared  with  similar  Dirac-Slater  and  Hart ree- Fock  calculations,  and  with  experi¬ 
ment.  It  is  found  that  the  calculated  binding  energies  are  generally  in  very 
good  agreement  with  experiment  and  are  superior  to  Dirac-Slater  and  Hartree-Foc 
results.  The  relativistic  contribution  is  shown  to  be  significant  even  for 


lighter  atoms. 
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Relativistic  calculations  of  atomic  core  binding  energies  are  now  being 
carried  out  with  increasing  accuracy.  It  has  been  established  that  relativistic 
effects  can  play  a  major  role  in  determining  the  magnitudes  of  these  energies. 
Comparisons  with  accurate  experimental  values  remain  sparse  however,  partly  be¬ 
cause  measurements  on  ..toms  have  until  recently  been  essentially  restricted  to 
the  rare  gases.  This  situation  is  now  changing  due  to  the  availability  of  data 
obtained  by  a  number  of  experimental  techniques.^ 

In  this  study  we  have  selected  one  subshell  in  potassium,  rubidium, 

cesium,  magnesium,  zinc  and  cadmium  and  computed  the  binding  energy  using  the 

2 

Dirac- Fock  program  of  Desclaux,  as  the  difference  between  the  total  energy  of 
the  hole  state  and  the  neutral  ground  state  ("ASCP'  method).  Where  appropriate 
we  have  calculated  explicitly  the  energies  of  electronic  states  of  the  ions  cor¬ 
responding  to  varying  .1  values,  although  for  these  cases  the  splitting  is  not 

observable  experimentally.  The  choice  of  systems  for  this  study  was  dictated  by 

3  4 

the  availability  of  experimental  values  from  photoelectron  spectroscopy.  '  For 
comparison,  we  have  also  calculated  the  binding  energies  of  the  corresponding 
levels  in  the  rare  gases,  although  most  of  these  binding  energies  have  been  com¬ 
puted  previously  by  other  authors  to  a  similar  degree  of  accuracy.  Finally,  we 
also  report  the  corresponding  binding  energies  in  the  cations  of  the  metal  atoms. 
Our  calcul?tions  are  compared  with  other  theoretical  results  and  with  experiment 
in  an  attempt  to  determine  the  magnitude  of  the  relativistic  contribution  in 
each  case.  The  data  are  shown  in  Tables  1  and  2. 

3  4 

First  we  compare  our  ASCF  results  with  experiment.  ’  The  agreement  is 
very  good  in  general,  with  most  of  the 

discrepancies  being  0.7  eV  or  less.  Notable  exceptions  are  the  neighboring  ele¬ 
ments  Kr  and  Rb,  where  the  calculated  values  are  too  high  by  ~3eV. 
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This  is  certainly  not  snrpri  ing  for  the  krypton  3p  since  it  is 

well  known  that  the  s i ngle-pa rt i c le  description  breaks  down  in  this  case  due  to 

the  presence  of  close-lying  doubly- ionized  states  with  energies  close  to  the  3p  hole 

states.^  Furthermore,  although  these  effects  are  thought  to  cease  beyond  2  = 

36,  there  is  some  evidence^  that  the  rubidium  3])  region  may  exhibit  similar 
features.  For  the  majority  of  the  levels  reported  here,  the  Dirac-Fock  binding 
energy  is  lower  than  the  experimental  value.  This  is  to  be  expected  since  the 
primary  reason  for  the  discrepancy  is  the  difference  between  the  correlation 
energy  contribution  to  the  ground  (neutral)  state  and  the  hole  state  energies. 

In  most  cases,  more  correlation  is  expected  in  the  n-electron  system  than  in  the 
n-1  system;  this  results  in  an  increase  in  the  binding  energy  when  the  correla¬ 
tion  energy  is  included. 

A  number  of  approximate  relativistic  calculations  have  been  reported, 
the  most  extensive  being  the  work  of  Huang,  et  al.  using  the  Slater  approxi¬ 
mation  for  the  exchange  integrals.  The  method  used  in  the  present  work  differs 
from  that  of  Ref.  6  primarily  in  that  exchange  is  treated  exactly.  It  is  of  in¬ 
terest  therefore  to  compare  the  binding  energies  obtained  by  the  two  methods. 

This  is  done  in  columns  2  and  4  of  Table  1  and  columns  2  and  3  of  Table  2.  The 
Dirac-Slater  results  are  in  reasonable  agreement  with  experiment,  but  generally 
the  agreement  is  not  as  good  as  with  Dirac-Fock,  as  expected.  Interestingly,  the 
Dirac-Slater  binding  energies  are  invariably  lower  than  the  Dirac-Fock  binding 
energies.  This  is  perhaps  due  to  an  overestimation  of  the  exchange  contribution 
in  the  hole  state  compared  to  the  ground  (neutral)  state. 

Another  reason  for  carrying  out  these  calculations  is  to  compare  with 
analogous  nonrelativistic  results.  Recently  Broughton  and  Bagus  reported  ex¬ 
tensive  binding  energy  results  for  selected  levels  of  a  large  number  of  atoms 
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and  ions.  Neither  our  calculations  nor  theirs  suffer  from  any  basis-set  defi¬ 
ciencies,  since  numerical  methods  were  used  in  both  cases;  therefore  any  differ¬ 
ences  must  be  attributed  to  relativistic  effects.  A  dramatic  example  of  this 
is  the  Mg  Is  level  where  the  difference  amounts  to  M.9  eV.  The  accepted  expla¬ 
nation  is  that  relativity  causes  the  s  levels  to  contract  thus  raising  the  bind- 
8 

ing  energy.  All  the  other  levels  we  wish  to  compare  with  the  results  of  Ref.  7 
are  spin-orbit  split,  an  effect  that  is  of  course  absent  in  the  non-relativistic 
calculations.  The  single  level  obtained  in  Hartree-Fock  does  not  truly  corres¬ 
pond  to  an  average  of  the  spin-orbit  components,  as  can  be  seen  from  examination 
of  the  relativistic  hamiltonian.^  Using  the  non-relativistic  hamiltonian  as  a 
reference,  there  is  no  reason  to  expect  that  the  two  spin-orbit  levels  should 
be  affected  symmetrically  by  the  addition  of  all  the  extra  terms  which  account  for 
relativistic  effects.  We  should  therefore  compare  the  non-relativistic  binding 
energy  to  both  spin-orbit  split  levels  rather  than  their  mean.  For  p  levels,  we 
note  that  the  non-relativistic  result  is  always  closer  to  the  binding  energy  of 


the  P3/2  leve*  while  for  d  levels  it  more  closely  resembles  the  d^2  level.  The 
trend  in  d  levels  has  been  observed  previously  in  our  study  of  cesium  and  iodine, 
and  their  ions.*1'  Its  origin  is  simply  that  in  the  relativistic  case  the  s  and 


p  electrons  are  more  tightly  bound  than  in  the  non-relativistic  approximation, 
while  for  d  electrons  the  opposite  is  true.  Thus  the  less  bound  of  the  p  and 
more  tightly  bound  of  the  d  spin-orbit  conponents  lie  nearest  the  Hartree-Fock 
values.  The  relativistic  shift  of  the  spin-orbit  pair  is  often  about  half  the 
spin-orbit  splitting,  resulting  in  a  near  coincidence  of  the  Hartree-Fock  and  one 
of  the  Dirac -Foci,  values.  Note  also  that,  as  expected,  the  non-relativistic 
results  do  not  agree  as  well  with  experiment  as  do  the  relativistic  ones,  with 
discrepancies  of  2-3  eV  in  some  cases.  Particularly  interesting  is  the 
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nonrelativistic  value  for  the  Cs  3d,  which  is  about  4  eV  higher  than  the 
experimental  binding  energy.  Similar  behavior  is  expected  for  Xe  3d.  The 
Cd  3d  levels  show  the  same  effect  but  to  a  somewhat  lesser  extent;  here  a  dif¬ 
ference  of  about  2  eV  is  observed.  In  all  cases  the  orbitals  either  expand  or 
contract  in  comparison  to  their  non-relativistic  approximations  depending  on 
their  £  values,  as  noted  above.  Thus  care  should  be  exercised  in  explaining  the 
differences  between  non-relat ivistic  Hartree-Fock  binding  energies  and  experi¬ 
ment  since  our  results  clearly  show  that  the  correlation  energy  contribution  is 
not  the  only  source  of  disagreement  even  for  these  lighter  atoms.  We 
would  expect  this  to  be  the  case  in  molecules  containing  these  atoms. 

We  have  also  calculated  core  binding  energies  of  the  corresponding 
positive  ions,  primarily  for  comparison  wi  the  results  of  Broughton  and  Bagus.7 
Although  the  relativistic  and  non-rel ativi c  atom-ion  shifts  are  close,  they 
are  really  not  quantitatively  the  same,  indicating  that  relativistic  effects  do 
not  cancel  out  completely  in  going  from  the  neutral  to  the  ion.  In  general,  the 
relativistically  calculated  shift  is  high-  r,  sometimes  by  as  much  as  0.2  eV. 

This  is  most  likely  due  to  the  fact  that  the  valence  electrons,  being  of  s  sym¬ 
metry,  are  more  contracted  in  the  relativistic  case  and  therefore  more  effective 
in  screening  core  electrons  from  the  nucleus. 

Turning  to  the  relaxation  energies  (the  difference  between  the  negative 
of  the  orbital  energy  and  the  ASCF  energy)  we  note  that  spin-orbit  components 
possess  different  relaxation  energies.  This  has  been  noted  once  before10  as 
being  clearly  due  to  the  added  relaxation  contribution  of  the  1  ♦  4  component  to 
ionization  from  the  1-4  level  but  not  the  reverse.  We  also  note  that  the  non¬ 
relativistic  relaxation  energies  are  always  lower  than  the  relaxation  energies 
of  either  of  the  corresponding  spin-or*  r  components. 
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It  is  interesting  to  compare  the  relaxation  energies  in  isoelectronic 
species  such  as  the  alkali  metal  ions  and  the  rare  gases.  (The  results  for  the 
latter  are  shown  in  Table  2.)  In  every  case  the  lower-2  species  show  more  re¬ 
laxation.  Similar  behavior  was  noted  for  the  isoelectronic  ions  Cs*  and  I".*** 

The  explanation  probably  lies  in  the  more  diffuse  nature  of  the  orbitals  of  the 
lower-Z  species  thus  allowing  for  more  relaxation.  An  interesting  example  of 
this  is  Mg*"  and  Ne,  where  the  relaxation  energy  is  about  1  eV  less  in  Mg^*  than 
in  Ne. 

The  calculated  core  binding  energies  for  the  rare  gases  are  in  good 
agreement  with  experiment  except  for  krypton  as  noted  above.  Results  very  close 
to  ours  have  been  obtained  by  other  workers.** 

In  conclusion,  we  have  shown  that  single  configuration  relativistic 
calculations  yield  binding  energies  in  very  good  agreement  with  experiment  for 
most  of  the  elements  considered  in  this  work  and  that  relativistic  contribution 
can  be  significant  even  in  light  atoms. 
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TABLE  1 

Core  Binding  Energies  and  Relaxation 


Level 

Binding  Energy 

D-Sb 

H-FC 

D-Fd 

K  2p1/2 

K  2p3/2 

303.00 

300 . 20 

301 . 32 

303.07 
303. 11 
300. 38 
300. 31 

K>  2p3/2 

308.43 

307.47 

Rb  3p1/2 

Rb  3P3/2 

257. 10 

247.86 

248.94 

257.22 

2S7.25 

248.07 

248.01 

5p3/2 

255. 36 

254.54 

Cs  3d3/2 

CS  3d5/2 

745.80 

731 . 88 

749.67 

745.88 

745.89 
732.08 
732.07 

Cs  3d5/2 

7SS.32 

737.88 

Mg  Is 

1310.58 

1309.24 

1311.10 

Mg2*  Is 

1332.36 

1334.25 

Zn  2Pl/2 

1051.43 

1031.60 

1051.73 

Zn  2p3/2 

1028.02 

1028.47 

Zn  2p3/2 

1054.10 

1051.18 

Cd  3d3/2 

417.88 

420.55 

418.07 

Cd  3d$^2 

411.12 

411.38 

Cd2*  3d5/2 

440.14 

431.48 

Energies  (eV) 


Ion -Atom 
Shi  ft f 

Relaxat ion 

ft 

Energy 

Expe 

H-F 

D-F 

303.2 

12.14 

12.52 

300.5 

7.13 

12.30 

“ 

(7.10) 

10.99 

11.15 

254.3 

9.24 

9.93 

245.4 

6.51 

9.53 

” 

(6.42) 

8.29 

8.56 

745. 6g 

18.42 

19.46 

731. 6g 

S .  81 

18.74 

- 

(5.65) 

17.61 

17.93 

1311.3* 

23.15 

25.02 

25.71 

“ 

(23.12) 

21.95 

22.65 

1052.1* 

29.66 

1028.9* 

22.71 

27.62 

28.73 

- 

(22.49) 

25.36 

26.49 

418.8 

16.80 

17.34 

412.0 

20.10 

16.89 

- 

(19.59) 

15.09 

15.20 
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TABLE  1  (continued) 

aThe  negative  of  the  orbital  eigenvalue  minus  ASCF  energy;  the  H-F  results  are 
from  Ref.  7.  The  L)-l  results  were  obtained  using  the  ASCF  energy  averaged  over 
multiplets . 

^Dirac-Slater  ASCF  results  from  Ref.  6. 

cHart ree- Fock  ASCF  results  from  Ref.  7. 

dSingle  configuration  ASCF  Dirac-Fock  calculations  using  the  program  in  Ref.  2. 
Where  two  energies  are  listed,  these  correspond  to  states  of  least  (upper)  and 
greatest  (lower)  J  quantum  number. 

eFrom  Ref.  3  for  Zn  and  Cd  and  Ref.  4  for  the  rest. 

fThe  binding  energies  used  to  compute  this  shift  are  averaged  over  multiplets 
for  the  atoms.  The  numbers  in  parentheses  are  non-relat ivistic  shifts  from 
Ref.  7. 

*The  experimental  values  from  Refs  3  and  4  were  reduced  by  0.2  eV  in  view  of 
the  most  recent  experimental  Nels  value  of  870.21  eV  (ste  Table  2);  this  is 
the  level  used  to  calibrate  the  spectra. 
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TABLE  2 

Core  Binding  Energies  and  Relaxation  Energies  for  the  Rare  Gases  (eV) 


Level 

Binding  Energy 

Relaxation  Energy 

D-S3 

l)-Fb 

Exp 

D-F 

Ne  Is 

869. IS 

869.47 

870. 21C 

23.55 

Ar  2p1/2 

249.95 

250.28 

250. 78d 

11.82 

Ar  2p3/2 

247.80 

248. 16 

248.63d 

11.63 

Kr  3p1/2 

224.82 

225.01 

221.8® 

9.55 

Kr  3p3/2 

216.81 

217.03 

214. 2e 

9.17 

Xe  3d3/2 

688.81 

689.03 

689. 3Se 

19.11 

Xe  3d5/2 

676.16 

676.47 

676.70® 

18.44 

aFrom  Ref.  6. 

^Present  work, 
c  ® 

From  H.  Agren,  J.  Nordgren,  L.  Selander,  C.  Nordling,  and  K.  Siegbahn, 

J.  Electron  Spectrosc.  14^  27  (1978). 
d  0 

From  J.  Nordgren,  H.  Agren,  C.  Nordling,  and  K.  Siegbahn,  Phys.  Scr. 
19,  5  (1979). 

®From  Ref.  5. 
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Symmetries  of  Transition  States  in  Isoenergic  Reactions 


Sir: 

Several  authors  have  shown  that  transition  states  (ts)  of  chemical  reactions  must 
possess  unique  symmetry  properites.  A  systematic  set  of  rules  have  been  put  forth 
by  Mclver  and  Stanton1,  who  stated  for  example  that  a  transition  vector  can  not 
correspond  to  a  degenerate  vibrational  mode.  However  the  special  case  of  isoenergetic 
reactions,  in  which  products  and  reactants  have  the  same  energy  to  a  high  level  of 
approximation,  gives  rise  to  an  additional  and  powerful  set  of  symmetry  constraints 
which  have  apparently  not  previously  been  considered.  Most  reactions  are  not,  of 
course,  isoenergetic.  Perhaps  the  most  common  examples  of  this  type  of  process 
are  isotopic  exchange  reactions  and  molecular  rearrangements  between  energetically 
equivalent  conformations.  In  this  paper  the  symmetry  condition  that  governs  isoener- 
getic  reactions  is  derived  and  applied  to  three-  and  four-  center  mechanisms.  Specific 
application  is  then  made  to  representative  examples  of  the  two  kinds  of  reaction 
mentioned  above. 

It  will  be  assumed  throughout  that  the  transition  state  corresponds  to  a  maxi¬ 
mum  in  the  energy  along  the  reaction  coordinate,  and  is  unique.  It  is  then  straightfor¬ 
ward  to  show  that  the  point  group  of  the  nuclei  in  the  TS  must  contain  at  least  one 
symmetry  operation.  If  products  and  reactants  are  of  equal  energy,  the  forward  and 
reverse  directions  along  the  reaction  coordinate  must  be  energetically  equivalent.  The 
vibrational  mode  of  the  TS  that  represents  reaction  in  these  two  directions  must  be 
symmetric,  so  that  the  Taylor  series  expansion  of  the  energy  in  terms  of  the  reaction 
coordinate  may  have  no  odd  contributions.  Odd  terms  may  appear  only  if  they  are 
totally  symmetric  under  all  the  symmetry  operations  of  the  group2.  However  if  the 
TS  had  no  symmetry,  then  every  vibration  would  be  totally  symmetric.  Hence  the 
vanishing  of  odd  terms  in  the  energy  implies  at  least  one  symmetry  operation. 

Not  every  vibration  of  the  TS,  regardless  of  symmetry,  will  necessarily  lie  along 
the  reaction  coordinate.  It  is  necessary  to  also  specify  what  groups  of  atoms  represent 
products  and  reactants.  In  the  case  of  isoenergetic  reactions  it  is  clearly  necessary  that 
there  be  two  ways  of  grouping  atoms  as  products  and  reactants  that  are  symmetrically 
equivalent,  so  that  the  vibrations  separating  them  will  be  energetically  equivalent. 
For  example  consider  the  reaction  A1-A2  +  A3  — ♦  Ai  +  A2-A3.  Symmetrically 
equivalent  groupings  in  the  TS  are  [(Ai,A2),(As)]  and  |(Ai),(A2,A3)].  Thus,  the  bonds 
A1-A2  and  A2-A3  must  be  equal  and  the  TS  has  C2v  symmetry.  This  is  also  what 
one  would  expect  intuitively. 

Now  consider  a  reaction  involving  four  centers.  In  this  case  the  result  is  less  intui¬ 
tively  apparant.  For  the  reaction  A1-A3  +  A2-A4  — »  A1-A2  +  A3-A4,  symmetrically- 
equivant  groups  of  centers  Ai,  A2,  A3  and  A4  are  [(Ai,A2),(A3,A4)],  |(Ai)A3),(A2,A4)j 
and  ((Ai,A4),(A2,A3)].  A  minimum  of  two  of  these  groups  must  be  symmetrically 
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equivalent,  say  the  first  two.  In  this  case  bond  A1-A2  and  A1-A3  are  equal,  as  are 
A3-A4  and  A2-A4  (or  equivalently  A1-A2  and  A2-A4,  A3-A4  and  Ai-A.s).  Again,  a 
TS  of  C2v  symmetry  is  obtained,  with  the  vertical  symmetry  plane  passing  through 
precisely  two  atoms.  It  is  important  to  note  however  that  the  symmetry  of  this  i  s 
is  not  unique  to  this  reaction.  Thus  for  the  reaction  A1-A2-A3  +  A4  — *  Aj  + 
A2-A3-A4  there  are  four  possible  symmetrically  equivalent  groups:  [(Ai,A2,A3),A4], 
[(Ai,A2,A4),A3j,  etc..  Two  of  these  groups  must  be  equivalent,  say  the  two  given. 
Now  there  are  two  distinct  possibilities:  A1-A2  =  A2-A4  and  A2-A3  =  A1-A2;  or 
A2-A3  =  A2-A4  and  A1-A2  *  A1-A2.  It  may  be  seen  that,  depending  on  the  way 
in  which  the  equivalent  groups  are  chosen,  a  IS  of  either  C%v  or  C2  symmetry  may 
result. 

An  example  of  an  isoenergetic  reaction  that  has  been  widely  studied  is  H2  +  D2 
-♦  2HD  which,  on  the  basis  of  experimental  data3,  has  been  thought  to  proceed  via  a 
bimolecular  mechanism.  The  symmetries  of  the  possible  paths  for  this  reaction  have 
been  studied  by  Wright4  and  Silver5,  who  however  did  not  consider  its  isoenergetic 
properties.  A  large  number  of  putative  TS  constrained  to  be  of  a  particular  symmetry 
have  been  computed  theoretically,  most  of  which  possess  the  C2v  symmetry  required 
by  the  above  arguments.  These  include  the  square6’7,8,  tetrahedron6’7’10,  rhombus6, 
kite6'8,  Ys  and  T9.  However  the  rule  derived  above  excludes  two  of  the  conformations 
which  are  the  most  favorable  energetically:  the  trapeaoid6'8  and  linear  chain6’'’11. 
Also  excluded  are  the  rhomboid6,7  and  rectangle8  among  others,  and  the  completely 
asymmetric  case5.  Another  interesting  example  is  the  inter  conversion  between  the 
two  poBBibile  rectangular  conformations  of  cyclobutadiene  and  its  singly-charged  ions. 
This  process  is  currently  thought  to  take  place  in  the  plane  of  the  molecule  and  to 
pass  through  a  square  conformation^  although  few  other  alternatives  been  studied 
thuB  far.  For  rearrangements  of  this  type  the  same  rule  applies  and  the  IS  must  have 
at  least  kite  symmetry. 

The  requirement  that  there  be  symmetrically  equivalent  groups  in  the  TS  of  an 
isoenergetic  reaction,  along  with  other  rules  governing  nuclear  and  orbital  symmetry, 
are  fairly  restrictive  and  may  uniquely  determine  the  TS  in  many  cases.  In  comput¬ 
ing  activation  energies  by  a b  initio  methods,  one  need  only  optimise  a  molecular 
geometry  of  the  minimum  symmetry.  If  a  bound  species  is  obtained  in  this  way,  it 
provides  a  direct  estimate  of  the  activation  energy.  If  a  bound  species  is  not  obtained, 
it  implies  that  the  proposed  mechanism  is  incorrect.  . 

Carl  S.  Ewig 

Department  of  Chemistry,  Vanderbilt  University 

Nashville,  Tennessee  37204 
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Abstract 


A  technique  for  determining  the  symmetries  of  transition  states  of  isoenergetic 
reactions  is  proposed,  and  applied  to  the  possible  three-  and  four-center  mechanisms. 
Specific  examples  given  are  the  exchange  of  a  H  atom  between  two  Hj  molecules,  and 
the  rearrangement  of  the  cyclobutadiene  molecule. 
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